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ABSTRACT 



> 

in 

*vj ' We present a study of the relation between dark matter halo mass and the baryonic content 

T-j- I of their host galaxies, quantified through galaxy luminosity and stellar mass. Our investiga- 

tion uses 154 deg^ of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lens- 
ing and photometric data, obtained from the CFHT Legacy Survey. To interpret the weak 
lensing signal around our galaxies we employ a galaxy-galaxy lensing halo model which 
allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to 
lenses at redshifts between 0.2 and 0.4, split into a red and a blue sample. We express the 
relationship between dark matter halo mass and baryonic observable as a power law with 

k> ' pivot points of 10^^ hjQ Lq and 2 x 10^^ h^Q Mq for luminosity and stellar mass respec- 

%^ . tively. For the luminosity-halo mass relation we find a slope of 1.56^gQg and a normaUsation 

. . .' of 1.26^006 ^ 10^"^ ^70^ -^0 for red galaxies, while for blue galaxies the best-fit slope is 

0.73to.o8 and the normahsation is 0.16 ± 0.03 x lO^^ /i^ ^^e- Similarly, we find a best-fit 
slope of 1.49to;o4 and a normalisation of 1.30lo;o9 ^ 10^^ Ko -^o f°r "^^e stellar mass-halo 
mass relation of red galaxies, while for blue galaxies the corresponding values are 0.831q q4 

and 0.88 ±0.13 x lO^'^ h^^^ Mq. All numbers convey the 68% confidence limit. For red 
lenses, the fraction which are satellites inside a larger halo tends to decrease with luminosity 
and stellar mass, with the sample being purely satellites for a stellar mass of 2 x 10*^ h^Q Mq. 
The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or 
stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, 
is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered 
environments than red galaxies. We also find that the halo model, while matching the lensing 
signal around red lenses well, is prone to over-predicting the large-scale signal for faint and 
less massive blue lenses. This could be a further indication that these galaxies tend to be more 
isolated than assumed. 

Key words: cosmology: observations - gravitational lensing: weak - galaxies: haloes - dark 
matter 
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1 INTRODUCTION 

In order to fully understand the mechanisms behind galaxy forma- 
tion, the connection between galaxies and the extensive dark mat- 
ter haloes in which they are enveloped must be studied in exhaus- 
tive detail. In pursuit of this precision, reliable mass estimates of 
both the baryonic and the dark matter content of galaxies are re- 
quired. The visible component may be evaluated using galaxy prop- 
erties such as the luminosity or the stellar mass, propertie s which 
can be derived via stellar synthesis m odels (Kauffmann et al. 2003; 
iGallazzi et al.ll2005l : lBell & de Jongil2001. : .Salim et al...2007.) . The 
dark matter, on the other hand, cannot be observed directly but 
must be examined through its gravitational influence on the sur- 
roundings. At the largest scales reached by haloes, optical tracers 
such as satellite galaxies are scarce. Furthermore, estim ates of halo 
mass from satellite galaxy kinematics (see, for example. iMore et al.l 
I2OI 11) do not only require spectroscopic measurements of a very 
large number of objects, which are unfeasible with current instru- 
mentation, but they also require the application of the virial theo- 
rem with all its associated assumptions. To study any and all galax- 
ies it is therefore desirable to use probes independent of these trac- 
ers, and independent of the physical state of the halo, but with the 
power to explore a large range of scales. These requirements are all 
satisfied by weak gravitational lensing. 

Gravitational lensing is a fundamental consequence of grav- 
ity. As light from distant objects travels through the Universe it 
is deflected by intervening matter. This deflection causes the dis- 
tant objects, or sources, to appear distorted (and magnified). In 
the weak regime the distortion is minute, and only by studying 
the shapes of a large number of sources can information about the 
foreground gravitational field be extracted. By examining the av- 
erage lensing distortion as a function of distance from foreground 
galaxies, or lenses, the density profiles of their dark matter haloes 
may be directly investigated; this technique is known as galaxy- 
galaxy lensing. First detected by iBrainerd et al.l d 1 996b . the field 
of galaxy-galaxy lensing has been growing rapidly, with increas- 
ing precision as survey area grows. Our understanding of the un- 
derlying physics also increases as the interpretation of the signal 
becomes more sophisticated. Simulations predict that dark mat- 
ter hal oes are well approximated by N avarro-Frenk-White profiles 
(NFW; Navarr o. Frenk. & Whitell 1996b and comparing such a pro- 
file to the observed signal around isolated lenses results in halo 
mass estimates. Galaxies and their haloes are not generally isolated, 
however, but reside in clustered environments. The ramification is 
that the interpretation of the observed galaxy-galaxy lensing sig- 
nal around foreground lenses becomes more complicated since the 
signal from neighbouring haloes also influences the result. To ad- 
dress this problem a number of approaches have been employed. 
Early studies modelled the lensing signal by associating all matter 
with galaxies and comparing the resulting shear field to the obser- 
vations in a rnaximum-Iikelihood approa ch JSchneider & Rixl 19971 ; 
[Hudson et alJl998l ; lHoekstra et al.l2004h . In this case the clustering 
of galaxies was acco unted for through the observed positions and 
Hudson et alj ([1998) explicitly attempted to correct for the offset 
signal seen by satellite galaxies in larger haloes. It was, however, 
an approximate description. Alternatively th e issue can be circum - 
vented by selecting only isolated lenses (see lHoekstra et alj|2005[) . 
This inevitably leads to a large reduction in the number of lenses, 
and the sample is no longer representative as it does not probe the 
full range of environments. 

Over the past decade a new app roach has gaine d trac- 
tion: the weak lensing halo model l lCoorav & ShethI 120021 ; 



Guzik & SeliakI I2OO2I; iMaridelbaum et"^ l2005bl ; [van Uitert et al] 
201 It lLeauthau d et al.l201 II) . Within the halo model framework, all 



haloes are represented as distinct entities, each with a galaxy at the 
centre. Enclosed in each main halo are satellite galaxies surrounded 
by subhaloes. In this work we seek to employ the halo model to gain 
a more accurate picture of galaxy-size dark matter haloes, allowing 
for a more precise analysis of the link between galaxies and the 
dark matter haloes they reside in. For this purpose we use image 
data from the completed Canada-France-Hawaii Telescope Legacy 
Survey (CFHTLS), and weak lensing and photometric redshift cat- 
alogues produced by the Canada-Fran ce-Hawai i Telescope Lens- 
ing Survey (CFHTLen^lHevmans e t al. 2012; Miller et alJI 201 3l ; 
iHildebrandt et al.l 12012 1. This work improves on the preliminary 
galaxy-galaxy lensing analysis carried out using a small subset of 
the CFHTLS and a single-halo model fi t to the inner regions only 
dParker et al.l2007h . Furthermore, unlike ICoupon et al.l (120 12 ^ who 
studied the clustering signal of galaxies for the full CFHTLS-Wide 
to constrain the evolution in redshift of the stellar-to-halo mass re- 
lation, our analysis is based on galaxy-galaxy lensing, which can 
directly constrain the average halo mass of galaxies on small scales. 
Three recent studies to use the galaxy-galaxy lensing halo 
model to constrain these relations are Mandelbaum et al. (200a), 



Ivan Uitert et al.l (l20Ilh (hereaft er IVUllh and lUeauthaud et al.l 
( l2012n . lMandelbaum et alj J2006I) studied the halo masses of lenses 
from the full area of the f ourth data release of the Sloan Digi- 
tal Sky Survey (SDSS DR4; lAdelman-McCarthv et alj2006l) using 
a galaxy-galaxy lensing halo model. The SDSS is very wide, but 
also very shallow which means that for low luminosity galaxies it 
is highly powerful while it lacks the depth to constrain the halo 
masses of higher-luminosity galaxies which a re at hig her redshifts 



on average. A similar study was performed bv lVU III using an ear 



lier implementation of the halo model software used for this paper. 
That study exploited a 300 deg^ overlap between the SDSS DR7 
and the intermediate-depth second Red-sequence Cluster Survey 
(RCS2; Gilbank et al. 201 1). The SDSS data were used to identify 
the lenses, but the lensing analysis was performed on the RCS2, 
improving greatly at the high mass end on the previous analy- 
sis based on the shallow SDSS alone. However, while the I VU III 
lenses had accurate spectroscopic redshift estimates, there were no 
redshift estimates available for the sources at the time. Thus the 
work presented here has, aside from the increased depth down to 
i^B = 24.7, a further advantage over the IVUIll analysis owing 
to the high-precision phot ometric redshifts availab le for all objects 
used in our analysis (see JHildebrandt et al.ll20I2b . This makes it 
possible to cleanly separate lenses from sources and therefore min- 
imises the contamination by satellites. It also allows for optimal 
wei ghting of the lensing sig nal. 

iLeauthaud et al.l ( I2OI2I) combined several techniques to con- 
strain the relation between halo mass and stellar mass using data 
from the deep spac e-based Cosmic Evolution Survey (COSMOS; 
IScoville et alj|2007l) . They did not, however, refine their results by 
splitti ng their lens sample according to galaxy type. In a follow-up 
paper, [Tinker et all ( |2012|) did split the COSMOS sample into star 
forming and passive galaxies to study the redshift evolution of the 
same relation, but limited their study to massive galaxies located 
centrally in a group-sized halo. Thanks to the large area and depth 
of the CFHTLenS, we are in this paper able to investigate the rela- 
tion for blue and red galaxies separately without limiting our sam- 
ple in that way. We provide here a detailed comparison between 
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our results and those quote d in iMandelbaum et alj (|2006h, IVUllI 
and iLeauthaud et alj ( l2012h . but leave [Tinker et alJ ( 1201 2|) due to 
the large difference in sample selection between our analysis and 
theirs. 

This paper is organised as follows: we introduce the data 
in Section [2] and in Section [3] we review our halo model and 
the formalism behind it. We investigate the lensing signal as a 
function of luminosity in Section |4] and as a function of stel- 
lar mass in Section |5] In Section |6] we compare our results to 
the three previous studies introduced above, and we conclude 
in Section [7] The following cosmo logy is assumed through- 
out (WMAP7: fKomatsu et ai] l201lh : (Om, ^a, /i.crg, w) = 
(0.27,0.73,0.70,0.81,-1). 



2 DATA 

In this paper we present a galaxy-galaxy weak lensing analysis 
of the entire Wide part of the Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS-Wide). The unique combination of area 
and depth makes this survey ideal for weak lensing analyses. The 
CFHTLS was a joint project between Canada and France which 
commenced in 2003 and which is now completed. The survey area 
was imaged using the Megaprime wide field imager mounted at the 
prime focus of the Canada-France-Hawaii Telescope (CFHT) and 
equipped with the MegaCam camera. MegaCam comprises an ar- 
ray of 9 X 4 CCDs and has a field of view of I deg^. The wide 
synoptic survey covers an effective area of about 154 deg^ in five 
bands: u* , g', r', i' and z'. This area is composed of four inde- 
pendent fields, Wl^, each with an area of 23-64 deg^ and with a 
full multi-colour depth of j^b — 24.7 (source in the CFHTLenS 
catalogue). The images have been independently reduced within 
the CFHTLenS coll aboration, and for details on this data reduction 
process, we refer to lErben et ari(l2009ll2012l) . 

CFHTLenS has measured accurate shapes and photometric 
redshifts for 8.7 x 10'^ galaxi es JHevmans et alJl20I2l : lMiller et al] 
l20I3l : lHildebrandt et al.ll20I2h . The shear estimates for the sources 
us ed in this w o rk hav e been obtained using lensfit as detailed 



in iMiller et al.l ( 120131) , and thoroughly tested for system atics 
within the CFHTLenS collaboration (see .Hevmans et alj20I2h . All 
sources also have multi-band phot ometric redshift estimates as de- 
tailed in iHildeb randt et al. (2012). The catalogues we use in this 
wo rk are discussed in Heym ans et al. (2012), Mille r et al.l ( 120130 
and lHildebrandt et al.l ( l20I2h , with the exception of the stellar mass 
estimates. These estimates were obtained and tested for this paper 
and we therefore describe them in detail below. 



2.1 Stellar masses 

Our primary photometry analysis uses the Bayesian ph otomet- 
ric redshift software BPZ dSemtezllIoOOl : ICoe et alj|2006l) to esti- 
mate photometric redshifts after performing an extinction correc- 
tion on the mult i-colour m agnitudes. Using BPZ with a simple set 
of six modified IColeman et al. (1980) templates is our preferred 
met hod to estimate redshift s when using only five optical bands 
(see IHildebrandt et alj|2012n . and we note that it has been shown 
that the BPZ software is as accurate for photo metric redshift esti- 
mates as the alternative Bayesian LePHARq^ ( lAmouts et al.ll99^ : 
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Figure 1. CFHTLenS stellar masses compared to those from the CFHT 
WIRCam Deep Survey (WIRDS) as a function of i^g-magnitude (top) and 
redshift (bottom) for red (dark purple solid dots) and blue (light green open 
triangles) galaxies. The upper panels in each plot show the dispersion in log 
stellar mass and the lower panels show the bias of the CFHTLenS stellar 
mass estimates relative to the WIRDS stellar mass estimates. 



Illbert et al.ll2006t) software (see IHildebrandt et al.ll2010l for a com- 
parison). For physical parameters such as stellar mass estimates, 
however, our preferred method is to u se a more complex set of 
galaxy templates. Using LePhare and lBruzual & Chariot ( l2003h 
templates has been proven to be a robust method to estimate physi- 
cal parameters (see Ilbert et al. 2010) and so we choose to use Le- 
Phare to estimate stellar masses. For a consistent analysis we also 
compute rest-frame luminosities from the same spectral template 
as used for the stellar mass estimates. 

We derive our stellar mass estimates by fitting synthetic spec- 
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tral energy distribution (SED) templates while keeping the redshift 
fixed at the BPZ maximum likelihood estimate. The SED templates 
are based on the stellar population synthesis (SPS) package devel- 
oped by Bruzual & Chariot (2003) a ssuming a Chabrier (2003) ini- 
tial mass function (IMF). Following lllbert et al.l <2010f) . our initial 
set of templates includes 18 models using two different metallici- 
ties (Zi = 0.008 Z0 and Z2 = 0.02 Zq) and nine exponentially 
decreasing star formation rates oc e~''^, where t is time and r 
takes the values r = 0.1, 0.3, 1, 2, 3, 5, 10, 15, 30 Gyr. The fi- 
nal template set is then generated over 57 starburst ages ranging 
from 0.01 to 13.5 Gyr, and seven extin ction values ra nging from 
( 3.05 to 0.3 using a lCalzetti et^ J2000l) extinction law. lllbert et al.l 
( 120 id) investigated the possible sources of uncertainty and bias by 
comparing stellar mass estimates between methods. The expected 
differ ence between our e stimates and those base d on a Sa lpeter 
IMF ('Amouts et al."2007^, a "diet" Salpeter IMF dBelfcOOSh , or a 
Kroupa IMF iBorch et al. 2006) is -0.24 dex, -0.09 d ex, orOdex 
respec tively (see lllbert et al.l2010l) . In their Section 4.2. lllbert et al.l 
( 120101) further argue that the choice of extinction law may lead to 
a systematic difference of 0.14, and the choice of SPS model to 
a median difference of 0.13-0.15 dex, with differences reaching 
0.24 dex for massive galaxies with a high star formation rate. 

We determine the errors on our stellar mass estimates via the 
68% confidence limits of the SED fit, using the full probability dis- 
tribution function. However, since we fix the redshift these errors 
tell us only how good the model fit is, and do not account for un- 
ce rtainties in the photome tric redshift estimates (see Section 5.2 
of iHildebrandt et alj|2012n . To assess the stellar mass uncertainty 
due to photometric redshift errors we therefore compare our mass 
estimates to those of the CFHT WIRCam Deep Survey (WIRDS; 
iBielbv et alj|2012n . The WIRDS stellar masses were derived from 
the CFHTLS Deep fields with additional broad-band near-infrared 
data using the same method as described here. We are thus compar- 
ing our CFHTLenS stellar mass estimates to other estimates which 
are also based on photometric data, but which have deeper pho- 
tometry leading to a more robust stellar mass estimate. The addi- 
tional near-infrar ed data allows us to rely on these estimates up to 
a redshift of 1.5 dPozzetti et al.ll2007h . For our comparison we use 
a total of 134,290 galaxies in the overlap between the CFHTLenS 
and WIRDS data, splitting our sample into red and blue galaxies 
using their photometric type Tbpz- Tbpz is a number in the range 
of [1.0, 6.0] representing the best-fit SED and we define our red 
and blue samples as galaxies with Tbpz < 1.5 and 2.0 < Tbpz < 
4.0 respectively, where the latter captures most spiral galaxies. A 
colour-colour comparison confirms that these samples are well de- 
fined. In Figure[T]we show the comparison between our stellar mass 
estimates and those from WIRDS as a function of magnitude (top, 
with galaxies in the redshift range [0.2, 0.4]) and redshift (bottom, 
with galaxies in the magnitude range [17.0, 23.5]). 

For the range of lens redshifts used in this paper, 
0.2 ^ zicns ^ 0.4, the total dispersion compared to WIRDS is then 
~ 0.2 dex for both red and blue galaxies. The lower panel in the 
bottom plot of Figure [T] shows that for red galaxies our stellar 
masses are in general slightly lower than the WIRDS estimates, 
with the opposite being true for blue galaxies. For galaxies brighter 
than Jab ~ 18> both the dispersion and the bias increase d ue to 
biases in the redshift estimates (see Hildebra ndt et al.ll2012h . The 
bias and dispersion also increase rapidly at magnitudes fainter than 
*AB ~ 23, again due to redshift errors. 

We emphasise that this comparison with WIRDS quantifies 
only the statistical stellar mass uncertainty due to errors in the pho- 
tometric redshifts and due to our particular template choice. Since 
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Figure 2. Magnitude (left panel) and photometric redshift (right panel) dis- 
tributions of galaxies in the CFHTLenS catalogue. For the left panel we 
show all galaxies in the CFHTLenS, while for the right panel we limit our 
sample to magnitudes brighter than i^g = 24.7. The upper limit of lens 
(source) magnitude used is shown with a dark purple dotted (light green 
dashed) line in the left panel, while our lens (source) redshift selection is 
marked with dark purple dotted (light green dashed) lines in the right panel. 
Though the lens and source selections appear to overlap in redshift. sources 
ai'e always selected such that they are well sepai'ated from lenses in redshift 
(see Section|T2). Furthermore, close pairs are down-weighted as described 
in Section[3JT] 



the mass estimates from both datasets have been derived using iden- 
tical method and template set, the systematic errors affecting stellar 
mass estimates are not taken into account above. The uncertain- 
ties arising from the choice of models and dust extinction law adds 
0. 15 dex and 0. 14 dex respectively to the error budget, as mentioned 
above, resulting in a total uncertainty of ~ 0.3 dex. 



2.2 Lens and source sample 

The depth of the CFHTLS enables us to investigate lenses with a 
large range of lens properties and redshifts, which in turn grants us 
the opportunity to thoroughly study the evolution of galaxy- scale 
dark matter haloes. As discussed bv IHildebrandt et al.l ( 120121) . the 
use of photometric redshifts inevitably entails some bias in red- 
shift estimates, and also in derived quantities such as luminosity 
and stellar mass. Our analysis is sensitive even to a small bias 
since our lenses are selected to reside at relatively low redshifts 
of 0.2 ^ 2;ions ^ 0.4, where 2 is understood to be the peak of the 
photometric redshift probability density function, unless explicitly 
stated otherwise (see Figure |2j. Because our lensing signal is de- 
tected with high precision, we empirically correct for this bias us- 
ing the overlap with a spectroscopic sample as described in Ap- 
pendix lAll Throughout this paper, we then use the corrected red- 
shifts, luminosities and stellar masses for our lenses. For the full 
survey area we achieve a lens count of Anions = 1.1 x 10*'. 

We then split our lens sample in luminosity or stellar mass 
bins as described in Sections |4| and |5] to investigate the halo mass 
trends as a function of lens properties. Since we have access to 
multi-colour data, we are also able to further divide our lenses in 
each bin into a red and a blue sample using photometric type as de- 
scribed in Section l2n We also ensure that our lenses are brighter 
than i^B < 23 which corresponds to an 80% completeness of the 
spectroscopic redshift sample we use to quantify the redshift bias 
discussed above. The high completeness ensures that the spectro- 
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scopic sample is a good representation of our total galaxy sam- 
ple. The galaxy sample is dominated by blue late-type galaxies for 
i'ab > 22, however, and we are thus unable to perform a reliable 
redshift bias correction for red lenses at fainter magnitudes due to 
a lack of objects. We therefore exclude red lenses with i^B > 22 
while allowing blue lenses to magnitudes as faint as Iab ~ 23. 
This selection is also illustrated in Figure |2] 

To minimise any dilution of our len sing signal due to pho- 
tometric redshift uncertainties, we follow iLeauthaud et alj J2012r) 
and use only sources for which the redshift 68% confidence limit 
does not overlap with the lens redshift. We further ensure that the 
lens and source are separated by at least 0.1 in redshift space. 
To verify the effectiveness of this separation, we compare the 
source number counts ar ound our lenses to that around random 
points (as suggested by ISheldon et alj l2004l Section 4.1). This 
test shows no significant evidence of contamination. The source 
magnitude is only limite d by the maximum CFHTLenS analysis 
depth of ^AR ~ 24.7 (s ee iHevmans et aLll2012l : iMiller et al.ll2013l : 
iHildebrandt et al.ll2012h . Note that we do not apply a redshift bias 
correction to source redshifts as there is no existing spectroscopic 
redshift survey at these faint limits. While it is important to correct 
our lenses for such a bias since the derived baryonic observables 
such as luminosity and stellar mass depend strongly on redshift, it 
is less important for the sources as the lensing signal scales with 
the ratio Dis/Dg, where Ds and Dis are the angular diameter dis- 
tances to the source, and between the lens and source respectively. 
This ratio is insensitive to small biases in the source redshifts. Our 
source count for the full survey (excluding masked areas) is then 
NaouTce = 5.6 X 10®, corresponding to an effective source density 
of 10.6 arcmin"'^ wh ere we use the source density definition from 
iHevmans et al.l ( 12012 !. Equation 1). 

The high quality of the CFHTLenS shear measurements 
has been verified via a se ries of systematics tests presented in 
IHevmans et"al] ilOlA and iMiller et al.l J2013h . To further illus- 
trate the robustness of the shears we perform two separate analy- 
ses specifically designed to test the galaxy-galaxy lensing signal. 
First, we use a sample of magnitude-selected lenses across the en- 
tire surv ey and compare the resulting weak lensing signal to that 
found bv lParker et all l l2007h for a 22 deg^ subset of the CFHTLS 
data, and to that found bv lVUlll for RCS2. Both previous analy- 
ses use shear measure ment software based on the cla s s of m eth- 
ods first introduced by iKaiser. Squires. & Broadhursli ( 119950 and 
known as KSB. The details of the comparison may be found in Ap- 
pendix lBT] and we find that the signal we measure agrees well with 
these earlier shear measurements. The second test, as described in 
Appendix IB 2 1 uses the seeing of the images to test for any potential 
multiplicative bias still remaining. We find that this bias is consis- 
tent with zero. 



3 METHOD 

To analyse the dark matter haloes surrounding our lenses we use 
a method known as galaxy-galaxy lensing, and compare the mea- 
sured signal with a halo model. In this section we will introduce the 
basic formalism and give an overview of our halo model. 



3.1 Galaxy-galaxy lensing 

The first-order lensing distortion, shear, is a stretch tangentially 
about a lens, induced by the foreground structure on the shape 



of a background source galaxy. Assuming that sources are ran- 
domly oriented intrinsically, the net alignment caused by lens- 
ing can be measured statistically from large source samples. In a 
galaxy-galaxy lensing analysis, source galaxy distortions are aver- 
aged in concentric rings centred on lens galaxies. We measure the 
tangential shear, 7t, as a function of radial distance from the lens 
this way, and also the cross shear, 7x , which is a 45° rotated signal. 
When averaged azimuthally, the cross shear can never be induced 
by a single lens which means that it may be used as a systematics 
check. The amplitude of the tangential shear is directly related to 
the differential surface density AE(r) — E(< r) — E(r), i.e. the 
difference between the mean projected surface mass density en- 
closed by r and the surface density at r, via 



AE(r) = Ecrit(7tW> 

with Ecrit the critical surface density 



Ecrit — 



Ds 



AttGDiDi, 



(1) 



(2) 



where Di is the angular diameter distance to the lens, and Ds and 
Dis are defined as before. Here, c is the speed of light and G is the 
gravitational constant. By comparing differential surface densities 
rather than tangential shears, the geometric factor is neutralised and 
the amplitudes of the signals can be directly contrasted between 
different samples. The only caveat is that the properties of lenses 
depend on the lens redshift so this difference still has to be taken 
into account. 

We calculate the weighted average shear in each distance bin 
from the lens by combining the shear measurement weight w with 
the geometric lensing efficiency rj — {DiDis)/Ds as described in 
IVelander et al.l ( 1201 iL Appendix B.4). By using r; we down-weight 
close pairs and can minimise any influence of redshift inaccuracies 
on the measured signal that way. We quantify any remaining red- 
shift systematics by calculating a correction factor for each mass 
estimate based on the redshift error distribution; see Appendix IA2I 
for details on how this is done. The average shear, scaled to a ref- 
erence redshift, is then given by 



Mr)} 






(3) 



where 77/ — rj/rj-rd is the lensing efficiency weight factor with ?7rcf 
a reference lensing e fficiency value. Th e lensing weight Wi is de- 
fined in Equation 8 o fMiller et al.l (l2013i) . and accounts both for the 
ellipticity measurement error and for the intrinsic shape noise. Fi- 
nally, we convert the average shear to AE(r) using the Ecrit com- 
puted for the reference lens and source redshifts. 

The CFHTLenS shears are affected by a small but non- 
negligible multiplicative bias. Miller et alj ( 120131) have modelled 
this bias using a set of simulations specifically created as a 'clone' 
of the CFHTLenS, obtaining a calibration factor 7n(i^sN, ''gai) as 
a function of the signal-to-noise ratio, !/sn, and size of the source 
galaxy, rgai. Rather than dividing each galaxy shear by a factor 
(1 -\- m), which w ould lead to a biased calibration as discussed in 
IMiller et al.l ( l20I3h . we apply it to our average shear measurement 
in each distance bin using the correction 



The lensing signal is then calibrated as follows: 

(7(r)> 



(7°^'(r)> 



1 + K{r) 
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The effect of this correction term on our galaxy-galaxy analysis is 
to increase the average lensing signal amplitude by at most 6%. 
Though t here will be some uncertainty associated with this term, 
IKilbinger et alj (20130 find that it has a negligible effect on their 
shear covariance matrix. The calibratio n factor m enters linear ly in 
our Equation|5] while it is squared in the lKilbinger et alJ ( 120131) cor- 
relation function coiTection factor, thus amplifying its effect. The 
conclusion we draw is therefore that the impact of the calibration 
factor uncertainty will be insignificant in th is work. We also apply 
the additive c-term correction discussed in iHevmans et al.l ( 120120 
but find that it does not change our results either. 

The circular averaging over lens-source pairs makes this type 
of analysis robust against small-scale systematics introduced by for 
example PSF residuals in the shape measurement catalogues. Be- 
cause the galaxy-galaxy lensing signal is more resilient to system- 
atics than cosmic shear, we choose to maximise our signal-to-noise 
by using the full CFHTLenS area (except for masked areas) rather 
than removing the fields th at have not passed the cosmic shear sys- 
tematics test described in iHevmans et alj ( 120121) . However, there 
could be spurious large-scale signal present owing to areas being 
masked, or from lenses close to an edge, such that the circular av- 
erage does not cover all azimuthal angles. We correct for such spu- 
rious signal using a catalogue of random lens positions situated out- 
side any masked areas; the number of random lenses used is 50,000 
per square-degree field, which amounts to more than ten times as 
many as real lenses. The stacked lensing signal measured around 
these random lenses is evidence of incomplete circular averages 
and will be present in the observed stacked lensing signal as well. 
Because of our high sampling of this random points signal, we can 
correct the observed signal measured in each field by subtracting 
the signal around the random lenses. This random points test is dis- 
cussed in more detail in Mandelbaum et a l. ( 2005a). The test shows 
that for this data, individual fields do indeed display a signal around 
random lenses which is to be expected, even in the absence of any 
shapemeasurement error, due to cosmic shear and shot noise, and 
due to the masking effect mentioned above. Averaged over the en- 
tire CFHTLenS area the random lens signal is insignificant relative 
to the signal around true lenses ranging from ~ 0.5% to ~ 5% over 
the angular range used in this analysis. Additionally, to ascertain 
whether including the fields that fail the cosmic shear systematics 
test biases our results, we compare the tangential shear around all 
galaxies with 19.0 < iab < 22.0 in the fields that respectively 
pass and fail this test, and find no significant differences between 
the signals. 

3.2 The halo model 

To accurately model the weak lensing signal observed around 
galaxy-size haloes, we have to account for the fact that galaxies 
generally reside in clustered environments. In this work we do this 
by employing the halo model software first introduc ed inlVUlll . 
For full details on the exact implementation we refer to lVUllI ; here 

we give a qualitative overview. 

Our h alo model builds on work p resented in lOuzik & SeliakI 

( l2002h and lMandelbaum et al.] ( l2005bl) . where the full lensing sig- 
nal is modelled by accounting for the central galaxies and their 
satellites separately. We assume that a fraction ( 1 — a) of our galaxy 
sample reside at the centre of a dark matter halo, and the remaining 
objects are satellite galaxies surrounded by subhaloes which in turn 
reside inside a larger halo. In this context a is the satellite fraction 
of a given sample. 

The lensing signal induced by central galaxies consists of two 
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Figure 3. Illustration of the halo model used in this paper. Here we have 
used a halo mass of ^7200 = 10^^ /i^q Mq, a stellar mass of M« = 
5 X 10^" /lyQ Mq and a satellite fraction of a = 0.2. The lens redshift is 
•^icns = 0-5- Dark purple lines represent quantities tied to galaxies which 
are centrally located in their haloes while light green lines correspond to 
satellite quantities. The dark purple dash-dotted line is the baiyonic com- 
ponent, the light green dash-dotted line is the stripped satellite halo, dashed 
lines are the 1-halo components induced by the main dark matter halo and 
dotted lines are the 2-halo components originating from nearby haloes. 



components: the signal arising from the main dark matter halo (the 
I-halo term AE^'^) and the contribution from neighbouring haloes 
(the 2-halo term AE^*^). The two components simply add to give 
the lensing signal due to central galaxies: 



AEco 



AE„„,,t + AE„„„ 



(6) 



In our model we assume that all main d ark matter haloes are well 
represented by an NFW density profile (^Navarr o. Frenk. & White! 
[ 1996) with a m ass-concentration relationship as given by 
iDuffv et al.l (20081) . The halo model parameters resulting from an 
analysis such as ours (see, for example. Section |4j are not very 
sensitive to the exact halo concentration, however, as discussed in 
IVUll! . To compute the 2-halo term, we use the non-linear power 
spectrum from Smith et al. ( 2003) . We also assume that the de- 
pend ence of the galax y bias on mass follows the prescription 
from Sheth_et ak! ilOOV), incorporating the adjustments described 
in Tinker et al. (2005). Note that this mass-bias relation is empiri- 
cally calibrated on large numerical simulations, and does not dis- 
criminate between different galaxy types. The bias prescription we 
adopt for our halo model does not in clude non-linear effects. Fig- 
ure 1 from lMandelbaum et al.l ( 120121) shows that non-linear bias af- 
fects the galaxy-dark matter cross-correlation coefficient at the 2% 
level at a comoving separation of 4 ft^,,^ Mpc. The magnitude of the 
effect diminishes with increasing distance to 1% at 10 h^Q Mpc, 
and the influence on our 2-halo terms should be comparable. The 
affected regime, where the 1-halo and 2-halo terms overlap, is no- 
toriously difficult to model however. One major issue is that of halo 
exclusion which attempts to account for the way neighbouring dark 
matter haloes overlap. Since the effect of non-linearity is likely 
small compared to other modelling uncertainties on these scales, 
and since the affected range is beyond that used to determine halo 
masses in this paper (see Section|4j, we choose not to include non- 
linear biasing in our model. Finally, we note that the central term 
essentially assumes a delta function in halo mass as a function of 
a given observable since we do not integrate over the halo mass 
distribution. For a given luminosity bin, for example, the particu- 
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lar mass distribution within that bin therefore has to be accounted 
for. We do correct our measured halo mass for this in the follow- 
ing sections, assuming a log-normal distribution, and the correction 
method is described in Appendices IA2l and lA3l for the luminosity 
and stellar mass analysis respectively. 

We model satellite galaxies as residing in subhaloes whose 
spatial distribution follows the dark matter distribution of the main 
halo. The number density of satellites in a halo of a given mass is 
described by the halo occupation distribution (HOD) which is com- 
monly par ameterised through a power law of the form ( A^ ) — M"^. 
Following iMandelbaum et al.l (20051?), we set e = 1 for masses 
above a characteristic mass scale, defined to be three times the typi- 
cal halo mass of a set of lenses. For masses below this threshold, we 
use e = 2. In our model, the subhaloe s have been tidally stripped 
of dark matter in the outer regions. As IMandelbaum et al.l COOSbl) 
did, we adopt a truncated NFW profile, choosing a truncation radius 
of 0.4r2oo beyond which the lensing signal is proportional to r^^, 
where r is the physical distance from the lens. This choice results 
in about 50% of the subhalo dark matter being stripped, and we 
acquire a satellite term which supplies signal on small scales. Thus 
satellite galaxies add three further components to the total lensing 
signal: the contribution from the stripped subhalo (AE'^*'^'^), the 
satellite 1-halo term which is off-centre since the satellite galaxy is 
not at the centre of the main halo, and the 2-halo term from nearby 
haloes. Just as for the central galaxies, the three terms add to give 
the satellite lensing signal: 



AE,,t = AEf;''' + aeL^ + aeL^ 



(7) 



There is an additional contribution to the lensing signal, not 
yet considered in the above equations. This is the signal induced by 
the lens baryons ( AE^^^ This last term is a refinement of the halo 
model presented in I VU 111 necessary since weak lensing measures 
the tot al mass of a system and not just the dark matter mass. Fol- 
lowing |Leauthaud_et_alJ ( 120111) we model the baryonic component 
as a point source with a mass equal to the mean stellar mass of the 
lenses in the sample: 



AE^ 



(M*> 



(8) 



This term is fixed by the stellar mass of the lens, and we do not fit 
it. Note that we choose not to include the baryonic term for neigh- 
bouring haloes, but its contribution is negligible. 

Finally, to obtain the total lensing signal of a galaxy sample 
of which a fraction a are satellites we combine the baryon, central 
and satellite galaxy signals, applying the appropriate proportions: 



AE: 



AE'^" + (1 



a)AEccnt + aAEsat 



(9) 



All components of our halo model are illustrated in Figure |3] In 
this example the dark matter halo mass is Af200 ~ 10^^ K^^Mq, 
the stellar mass is M* = 5 x 10^" /ifo^ -^j^o, the satellite fraction 
is a = 0.2, the lens redshift is 2:icns = 0.5 and Dis/Ds — 0.5. On 
small scales the 1-halo components are prominent, while on large 
scales the 2-halo components dominate. 

The halo model is necessarily based on a number of as- 
sumptions (such as the spatial distribution of satellites, the HOD 
choice, the mass-concentration relation, the bias prescription and 
the adopted cosmology parameters), some of which may be overly 
stringent or inaccurate. Ideally, uncertainties should be assigned to 
these assumptions and included in the final fit parameter error bud- 
get. Alas, reliable eiTor estimates for these assumptions are not yet 
available in literature. To fully investigate the model assumptions 
would require the use of high-resolution hydrodynamical Simula- 
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Figure 4. r'-band absolute magnitude distribution in the CFHTLenS cata- 
logues for lenses with redshifts 0.2 ^ Zions ^ OA (black solid histogram). 
The distribution of red (blue) lenses is shown in dotted dark purple (dot- 
dashed light green). Our lens bins ai'e marked with vertical lines. 



Table 1. Details of the luminosity bins. (1) Absolute magnitude range; (2) 
Number of lenses; (3) Mean redshift; (4) Fraction of lenses that are blue. 



Sample 


M,,W 


"Icns*^' 


(^>(3) 


.fbluc(^' 


LI 


[-21.0,-20.0] 


91224 


0.32 


0.70 


L2 


[-21.5,-21.0] 


33633 


0.32 


0.45 


L3 


[-22.0,-21.5] 


23075 


0.32 


0.32 


L4 


[-22.5,-22.0] 


12603 


0.32 


0.20 


L5 


[-23.0,-22.5] 


5344 


0.32 


0.11 


L6 


[-23.5,-23.0] 


1704 


0.31 


0.05 


L7 


[-24.0,-23.5] 


344 


0.30 


0.03 


L8 


[-24.5,-24.0] 


76 


0.30 


0.09 



tions, and that is beyond the scope of this work. We therefore note 
that the results we obtain are specific to the adopted halo model, 
and any differences with alternative models should be accounted 
for when comparing our results to other work. This is likely to have 
more of an effect on the satellite fraction than on the measured 
halo masses though, since the above assumptions mainly affect the 
model on scales beyond the range of the central 1-halo term. Addi- 
tionally, the errors on the derived halo masses and satellite fractions 
in this work may be slightly underestimated due to the exclusion of 
modelling uncertainties. 



4 LUIMINOSITY TREND 

The luminosity of a galaxy is an easily obtainable indicator of its 
baryonic content. To investigate the relation between dark matter 
halo mass and galaxy mass we therefore split the lenses into 8 bins 
according to IVIegaCam absolute r'-band magnitudes as detailed 
in Table [T] and illustrated in Figure |4] The lens property averages 
quoted in this and forthcoming tables are pure averages and do not 
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Figure 5. Galaxy-galaxy lensing signal around lenses which have been split into luminosity bins according to Table[T] modelled using the halo model described 
in Section [X2l The dark purple (light green) dots represent the measured differential surface density, AE, of the red (blue) lenses, and the solid line is the 
best-fit halo model. Triangles represent negative points that are included unaltered in the model fitting procedure, but that have here been moved up to positive 
values as a reference. The dotted error bars are the unaltered error bars belonging to the negative points. The squares represent distance bins containing no 
objects. For a detailed decomposition into the halo model components, we refer to Appendix Icl 



include the lensing weights, unless explici tly spec ified. The choice 
of bin limits follows the lens selection in lVUlll This choice will 
allow us to directly compare our results to the results shown in 
IVUllI since the RCS2 data have been obtained using the same fil- 
ters and telescope. We also split each luminosity bin into red and 
blue subsamples as described in Section \2A\ and proceed to mea- 
sure the galaxy-galaxy lensing signal for each sample, with errors 
obtained via bootstrapping 10'' times over the full CFHTLenS area, 
where the number of bootstraps ensure convergence of the mean. 
We then fit the signal between 50 h^Q kpc and 2 h^g Mpc with 
our halo model using a x^ analysis. Only the halo mass M200 and 
the satellite fraction a are left as free parameters while we keep 
all other variables fixed. When fitting, we assume that the covari- 
ance matrix of the lensing measurements is diagonal. Off-diagonal 
elements are general l y pres ent due to cosmic variance and shape 
noise, but lChoi et alj ( 120121) find that for a lens sample at a redshift 
range similar to that of our lenses the covariance matrix is diago- 
nal up to ~ 1 Mpc, which coiTesponds well to the largest scale we 
include in our fits (this is also confirmed via visual inspection of 
our matrices). Furthermore, Figure 7.2 from the PhD thesis of Jens 
Rodige^j shows that the off-diagonal elements are comparatively 
small. Hence we do not expect that the off-diagonal elements in 
the x^ fit will have a significant impact on the best-fit parameters. 
The results are shown in Figure |5] for all luminosity bins and for 
each red and blue lens sample, with details of the fitted halo model 
parameters quoted in Table |2] The halo masses in this table have 



http://hss.ulb.uni-bonn.de/2009/1790/1790.htm 



been corrected for various contamination effects as detailed in Sec- 
tion l4. H and Appendix|A] Note that the number of blue lenses in the 
two highest-luminosity bins, L7 and L8, is too low to adequately 
constrain the halo mass. In the following sections, these two blue 
bins have therefore been removed from the analysis of blue lenses. 

As expected, the amplitude of the signal increases with lumi- 
nosity for both red and blue samples indicating an increased halo 
mass. In general, for identical luminosity selections blue galax- 
ies have less massive haloes than red galaxies do. For the red 
sample, lower luminosity bins display a slight bump at scales of 
~ 1 hjQ Mpc. This is due to the satellite 1-halo term becoming 
important and indicates that a significant fraction of the galaxies in 
those bins are in fact satellite galaxies inside a larger halo. On the 
other hand, brighter red galaxies are more likely to be located cen- 
trally in a halo. The blue galaxy halo models also display a bump 
for the lower luminosity bins, but this feature is at larger scales 
than the satellite 1-halo term. The signal breakdown shown in Fig- 
ure |C2| ( Appendix let reveals that this bump is due to the central 2- 
halo term arising from the contribution of nearby haloes. We note, 
however, that in these low-luminosity blue bins, the model overes- 
timates the signal at projected separations greater than ^Ih^^Mpc. 
This could be an indicator that our description of the galaxy bias, 
while accurate for red lenses, results in too high a bias for blue 
lenses. Alternatively, the discrepancy may suggest that the regime 
where the 1-halo term transitions into the 2-halo term is not ac- 
curately described due to inherent limitations of the halo model, 
such as non-linear galaxy biasing, halo exclusion representation 
and inaccuracies in the non-linear matter power spectrum (see Sec- 
tion [3]2j. To optimally model the regime in question, the handling 



CFHTLenS: Relation between galaxy DM haloes and baryons 9 



Table 2. Results from the halo model fit for the luminosity bins. (1) Mean luminosity for red lenses [10^'' h^^ IjqI 

i^'q Mq]; (4) Best-fit satellite fraction for red lenses; (5) Mean luminosity for 
7) Scatter-corrected best-fit halo mass for blue lenses [10^^ h^g Mq]; 
(8) Best-fit satellite fraction for blue lenses. The fitted parameters are quoted with their la errors. Note that the blue results from the L7 and L8 bins are not 
used for fitting the power law relation in Section l4~n 



[10^*^ /lyg Mq]; (3) Scatter-corrected best-fit halo mass for red lenses [10 
blue lenses [10^" /i^q ^0]; (6) Mean stellar mass for blue lenses [10^ 



"ft; 
2 



' Lq]; (2) Mean stellar mass for red lenses 



Sample {Ll'='^)W (Mf'i)(2) 



^rcd(3) 



j,red(4) (LJ?'"°)(5) (Af^'"'=)(^) 



j^,^bluc(7) 



bluc{8) 



LI 
L2 
L3 

L4 
L5 
L6 

L7 
L8 



0.91 
1.74 
2.73 
4.28 
6.69 
10.4 
16.4 
25.4 



1.83 


3.08 


3.74 


10.3 


5.97 


15.5 


9.35 


34.3 


14.9 


60.7 


23.9 


123: 


35.6 


344; 


20.3 


762 



+0.98 

-0.95 

+2.37 

-1.73 

+2.72 

-2.31 

+6.00 

-4.43 

+5.90 

-6.64 

+ 15.2 

16.0 

+70.4 

89.6 

-320 

154 



1+0.05 
-0.04 

S5+004 
^■•30_o.04 
1+0.04 

-0.04 
1+0.04 

-0.05 
-1+0.00 

-0.03 
-1+0.00 

-0.04 
5+0.14 

-0.06 
-1+0.00 

-0.20 



0.54"^ 



0.31: 

0.21^ 
0.20J 
0.20J 
0.06^ 
O.20I 



1.08 

2.23 
3.52 
5.51 
8.44 
13.7 



0.50 


3.72 


1.10 


1.83 


1.83 


11.2 


3.00 


13.0 


4.63 


16.3 


7.88 


23.3 



+0.65 
-0.70 
+ 1.35 
-1.05 
+3.22 
-2.50 
+4.13 
-4.40 
+ 11.4 
-9.34 
+ 26.5 
-18.9 



0.00: 

0.07 
0.00 
0.00 
0.00: 
0.00 



0.02 

0.00 

+0.04 
0.04 

+0.03 
0.00 
0.06 
0.00 

+0.12 
0.00 

+0.14 
-0.00 
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Figure 6. Satellite fraction a and bias-corrected halo mass A/200 as a func- 
tion of r'-band luminosity. Dark purple (light green) dots represent the re- 
sults for red (blue) lens galaxies, and the dash-dotted lines show the power 
law scaling relations fit to the Figure|5]galaxy-galaxy lensing signal (rather 
than to the points shown) as described in the text. The dotted line in the 
lower panel shows the a prior applied to the highest-luminosity bins. 
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Figure 7. Constraints on the power law fits shown in Figure |6] In dark 
purple (light green) we show the constraints on the fit for red (blue) lenses, 
with lines representing the 67.8%, 95.4% and 99.7% confidence limits and 
stars representing the best-fit value. 



of these factors should perhaps be dependent on galaxy type, but 
that is not done here. The reason is that we do not currently have 
enough data available to investigate this regime in detail. In the fu- 
ture, however, it should be explored further. 



4.1 Luminosity scaling relations 

Before determining the relation between halo mass and luminosity 
we have to correct our raw halo mass estimates for two systematic 
effects. Firstly, we rely on photometric redshift estimates which do 
not benefit from the absolute accuracy of spectroscopic redshifts. 
We can therefore not be certain that a lens which is thought to be at 
a certain redshift is in fact at that redshift. If the redshift is different, 
then the derived luminosity will also be different which means that 



the lens may have been placed in the wrong bin. Though the lenses 
can scatter randomly according to their individual redshift errors, 
the net effect will be to scatter lenses from bins with higher abun- 
dances to those with lower abundances. The measured halo mass 
will therefore be biased. To correct for this effect we create mock 
lens catalogues and allow the objects to scatter according to their 
redshift error distributions. Secondly, the halo masses in a given lu- 
minosity bin will not be evenly distributed, which means that the 
measured halo mass does not necessarily correspond to the mean 
halo mass. The derivation of the factor we apply to our halo masses 
to correct for both these effects is detailed in Appendix I A2I 

The estimated halo masses for all luminosity bins, corrected 
for the above scatter effects, are shown as a function of luminosity 
in the top panel of Figure |6] Red lenses display a steeper relation- 
ship between halo mass and luminosity than blue lenses, and the 
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haloes of the blue galaxie s are in general less massive for a given 
luminosity bin. Following lVUl iL we fit a power law of the form 



M200 = Mo,. 



(10) 



with Lfid — 10 h^Q L^i Q a scaling factor chosen to be the r'- 
band luminosity of a fiducial galaxy. Rather than fitting to the final 
mass estimates we fit this relation directly to the lensing signals 
themselves (taking the scatter correction into account). We do this 
because the error bars are asymmetric in the former case, but the 
difference in results between the two fitting techniques is small. 
For our red lenses we find 



Mq.l = 1.261?,°^ X 10^=* /ifo^ Mq and /?l = lM+_°o°ot 



while for our blue lenses 
Mo.L =0.16 ±0.03 X 10^^ 



the corresponding numbers are 



hjo Mq 



and /3l = 0.73 



+0.09 



The 



parameters are quoted with their lo errors and the constraints for 
these fits are shown in Figure |7] Here we again see that the red 
lenses are better constrained than the blue. This is partly because 
we have more red lenses in most bins, and partly because red 
lenses in general are more massive at a given luminosity. 

The mass-to-light ratios, M2oo/{Lr), of our red sample 
range from 34tjJ /170 Mq LZ^, at the lowest luminosity bin to 
91I10 hro Mq Lq^ for L5. For our blue sample the numbers are 
34 ± 11 hro Mq L~^ for LI and 19 ± 2 hro Mq L~^ for L5. Be- 
yond L5, the mass-to-light ratio continues to increase, reaching 
300lgi® hro Mq Lq^ for red lenses in bin L8. In these highest lu- 
minosity bins a significant fraction of the red lenses maybe associ- 
ated with groups or small clusters, as pointed out bv lVUlll 



4.2 Satellite fraction 

The lower panel of Figure|6]shows the satellite fraction a as a func- 
tion of luminosity for both the red and the blue sample. At lower 
luminosities the satellite fraction is ~ 50% for red lenses and as 
the luminosity increases the satellite fraction decreases. This in- 
dicates that a large number of faint red lenses are satellites in- 
side a larger dark matter halo, consistent with previous findings 
(seelMandelbaum et al]|2006l : Ivan Uitert et alj|20 111 : ICoupon et d] 
I2OI2I) . In the highest luminosity bins the satellite fraction is difficult 
to constrain due to the shape of the halo model satellite terms (light 
green lines in Figure |3) becoming indistinguishable from the cen- 
tral 1-halo term (dark purple dashed), as discussed in Appe ndix Icl 
To ensure that our halo masses are not biased low we follow I VUl ll 
and apply a uniform satellite fraction prior to these bins, allow- 
ing a maximum a of 20%. This prior is marked in Figure |6] For 
blue lenses, the satellite fraction remains low across all luminosi- 
ties indicating that almost none of our blue galaxies are satellites, 
again consistent with previous findings. This may be a sign that 
blue galaxies in our analysis are in general more isolated than red 
ones for a given luminosity, a theory corroborated by the low sig- 
nal on large scales for blue galaxies (see Figure |C2] in Appendix Icl). 
Here we have made no distinction between field galaxies and galax- 
ies residing in a d enser environment; for a more in-depth study of 
this distinction see lGillis et al.l ( 120131) . 



5 STELLAR MASS TREND 

The galaxy luminosity as a tracer of baryonic content depends both 
on age and on star formation history. A galaxy's stellar mass does 
not have such dependence and may therefore be a better indicator 
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Figure 8. Stellar mass distribution in the CFHTLenS catalogues for lenses 
with redshifts 0.2 S^ ^lons ^ 0.4 (black solid histogram). The distribution 
of red (blue) lenses is shown in dotted dark purple (dot-dashed light green). 
Our lens bins are marked with vertical lines. 



Table 3. Details of the stellar mass bins. (1) Stellar mass range [/i^q Mq]; 
(2) Number of lenses; (3) Mean redshift; (4) Fraction of lenses that are blue. 



Sample 


logioM.*'' 


nions*^' 


(.)(3) 


/blue'*' 


SI 


[9.00,9.50] 


126406 


0.33 


0.981 


S2 


[9.50,10.00] 


78283 


0.32 


0.828 


S3 


[10.00,10.50] 


48957 


0.32 


0.391 


S4 


[10.50,11.00] 


37365 


0.32 


0.043 


S5 


[11.00,11.25] 


7474 


0.32 


0.003 


S6 


[11.25,11.50] 


2447 


0.31 


0.001 


S7 


[11.50,11.75] 


396 


0.30 


0.000 


S8 


[11.75,12.00] 


12 


0.31 


0.000 



of its baryonic content. In this section we study the relation be- 
tween galaxy stellar mass and the dark matter halo mass, dividing 
the lenses into 9 stellar mass bins as illustrated in Figure|8]with de- 
tails in Table [5] As we did for the luminosity analysis (Section lU 
we further split each stellar mass bin into a red and a blue sample 
using their photometric types to approximate early- and late-type 
galaxies. 

We measure the galaxy-galaxy lensing signal for each sample 
as before, and fit on scales between 50 hjg kpc and 2 h^Q Mpc 
using our halo model with the halo mass Af200 and the satellite 
fraction a as free parameters. Similarly to the previous section, the 
results are shown in Figure |9]for all stellar mass bins and for each 
red and blue lens sample, with details of the fitted halo model pa- 
rameters quoted in Table|4] There are no blue lenses available in the 
two highest stellar mass bins, and in bins S5 and S6 the number of 
blue lenses is too low to constrain the signal. We therefore remove 
them from our analysis in the following sections. 

The mean mass in each bin increases with increasing stellar 
mass as expected, resulting in an increased signal amplitude. Simi- 
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Figure 9. Galaxy-galaxy lensing signal ai'ound lenses which have been split into stellar mass bins according to Table |3] modelled using the halo model 
described in Section [X2l The dark purple (light green) dots represent the measured differential surface density, AS, of the red (blue) lenses, and the solid line 
is the best-fit halo model. Triangles represent negative points that are included unaltered in the model fitting procedure, but that have here been moved up to 
positive values as a reference. The dotted error bars are the unaltered error bars belonging to the negative points. The squares represent distance bins containing 
no objects. For a detailed decomposition into the halo model components, we refer to Appendix IdI 

Table 4. Results from the halo model fit for the stellar mass bins. (1) Mean luminosity for red lenses [10^*^ h^Q Lq]; (2) Mean stellar mass for red lenses 
[10^" h^Q Mq]; (3) Scatter-corrected best-fit mean halo mass for red lenses [10^^ h^g ^q]' (4) Best-fit satellite fraction for red lenses; (5) Mean lumi- 
nosity for blue lenses [10^" h^Q Lq]', (6) Mean stellar mass for blue lenses [lO^'' h^g Mq]; (7) Scatter-corrected best-fit mean halo mass for blue lenses 
[10^^ /lyQ Mq]; (8) Best-fit satelHte fraction for blue lenses. The fitted parameters are quoted with their la en'ors. Note that the red results from the SI and 
S2 bins, and the blue results from the S5 and S6 bins, are not used for fitting the power law relation in Section IsTTI 



Sample {1^.'"^)^ (Mf'l)'^) 



^red(3) 



a 



rcd(4) (LJ?'"°)(5) (M]?'"")^^ 



j\^bluo(7) 



„blue(8) 



SI 
S2 
S3 
S4 
S5 
S6 
S7 
S8 



0.22 
0.44 
1.06 
2.46 
5.38 
8.96 
14.3 
19.1 



0.23 


1.68 


0.66 


2.93 


1.97 


3.58 


5.65 


20.9 


13.0 


70.5 


22.6 


145: 


38.6 


459; 


62.7 


264 



+0.99 
-0.65 
+ 1.50 
-1.16 
+ 1.03 
-0.99 
+2.57 
-1.86 
+ 10.5 
-10.5 
+29.5 
-30.0 
+94.8 
-71.0 
+432 
-243 



1.00 
0.82 
0.51 
0.31 
0.30 
0.28 
0.20: 
0.20 



0.00 

0.03 
+0.07 

0.07 
+0.04 

0.04 
+0.03 

0.03 
+0.05 

0.06 
+0.09 

0.09 
+0.00 

0.07 
+0.00 

0.20 



0.41 
1.11 

2.87 
7.07 



0.18 


1 qfi+O-51 


0.54 


o qi+0.79 
'^•yj^-0.73 


1.59 


10 9+2-63 


4.29 


30. 9+ P.^ 



-+0.02 
-0.02 
00+°°2 

-+0.03 
-0.03 

-1+0.05 
-0.00 



0.05^ 



0.05: 

o.ooj 



lar to the luminosity samples in the previous section, the red lower- 
mass bins display a bump at scales of ~ 0.5 h~^ Mpc. Here the 
lowest bins contain less massive galaxies than the lowest luminos- 
ity bins and the bump is more pronounced, indicating that most of 
the galaxies in these low-mass samples are satellite galaxies. The 
contribution from nearby haloes is again clearly visible in the best- 
fit halo model of the lower-mass blue samples, though as noted in 
Section|4l this may be due to an inaccurate galaxy bias description 
for blue lenses. 



5.1 Stellar mass scaling relations 

Just as for the luminosity bins, we have to correct the halo mass 
estimates for two scatter effects: one due to errors in the stellar 
mass estimates and another due to halo masses not being evenly 
distributed within a given bin. We describe the correction for these 
effects in Appendix IA3I The best-fit halo masses, once corrected 
for these scatter effects, and satellite fractions a for each stellar 
mass bin are shown in Figure [TO] In the lower-mass bins, nearly 
all red lenses are satellites while for higher masses, the majority 
are located centrally in their halo. As discussed in Section l42l this 
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Figure 10. Satellite fraction a and halo mass AI200 as a function of stellar 
mass. Dark purple (light green) dots represent the results for red (blue) lens 
galaxies. Open circles show the points that have been excluded from the 
power law fit because of a high satellite fraction. The dotted line in the 
lower panel shows the a prior applied to the highest-stellar mass bins. 



<3a. 




Figure 11. Constraints on the power law fits shown in Figure 1101 In dark 
purple (light green) we show the constraints on the fit for red (blue) lenses, 
with lines representing the 67.8%, 95.4% and 99.7% confidence limits and 
stars representing the best-fit value. 



fraction is difficult to constrain for high masses due to the shape of 
the halo model satellite terms. We therefore apply the same uniform 
satellite fraction prior to the high-stellar mass bins as we did to the 
high-luminosity bins, allowing a maximum a of 20%. The overall 
low satellite fraction for blue galaxies, suggesting together with low 
large-scale signal that most blue galaxies are isolated, is consistent 
with the luminosity results. 



It is clear that the relation between dark matter halo and stellar 
mass is different for red and blue lenses. To quantify the difference, 
we fit a power law to the lensing signals in each bin simultaneously, 
similarly to our treatment of the luminosity bins in the previous 
section. The form of the power law is 

M200 = Mo,M ( -n— I (11) 



Mfi 



with Mfid ~ 2 X 10^^ h'^Q Mq a sc aling fa ctor chosen to be the 
stellar mass of a fiducial galaxy as in lVUlll . We note that for the 
lowest red stellar mass bins, though the halo model fits the data 
very well (see Figure|9), the sample consists of nearly 100% satel- 
lite galaxies as mentioned above. The central halo mass associated 
with these lenses is therefore effectively inferred from the satellite 
term, and thus constrained indirectly by the halo model and so we 
exclude the two lowest stellar mass bins from our analysis. 

The resuhing best-fit values for red lenses are 



1.301°°^ X 10^ 



rlMr: 



Ko IVIQ 



and Pm = 1.49 



+0.06 



Mo,M 

and for blue lenses Mo,m = 0.88 ± 0.13 x 10^^ /i^ Mq and 
Pm = 0.83lo'o4- W^ show the constraints and best-fit values in 
Figure [TT] The red lenses are clearly better constrained than the 
blue ones due to the stronger signal generated by these generally 
more massive galaxies. We note here that due to a lack of massive 
blue lenses in our analysis, the two galaxy type results probe 
different stellar mass ranges. The blue relation is limited to the 
low-stellar mass end only, while the red relation is constrained 
mostly at higher stellar masses. 

The baryon fraction, M«/M2oo, is fairly constant between 
stellar mass bins though it shows a tendency to decrease for red 



lenses from 0.055_ 



for S3 to 0.008lo;;52i, 



for S7. For blue 



lenses it conversely shows a slight increase from 0.009 



+0.005 



for 



^+0.002 



baryon conversion efficiency, though the particular environment 
each sample resides in affects the numbers. Since the red and blue 
samples probe different stellar mass ranges, we cannot directly 
compare the two. 



6 COMPARISON WITH PREVIOUS RESULTS 

Early galaxy-galaxy lensing based works that have investigated the 
relati on between lum i nosity and the virial mass of ga laxies in- 
clude [Ouzir&Seii^ ([2003) and lHoekstraetal.1 ilOok) . In these 
works, the mass is found to scale with luminosity as cxL^'* and 
j^^i.6±o.2^ respectively, in agreement with our findings. We focus, 
however, on comparing our halo mass results with those from three 
recent comprehensive galaxy-galaxy lensing halo model analyses 
which used data from thre e decidedly different surveys: the very 
wide but shallow SDSS ( M andelbaum et al . 2006), the moderately 
deep a nd wide RCS2 (VUllh and the very deep but narrow COS- 
MOS jLeauthaud et alj|2012l) . All four datasets are shown in Fig- 
ures [T2]and [13] with our results denoted by solid dots. 

We begin our comparison by noting that the various works 
emplo y diff erent g alaxy type separation criteria. lMandelbaum et al.l 
( 120061) and IVUll, base their selection on the brightness profile 
of the lenses, while we use the SED type. As both selection cri- 
teria are found to correlate well with the colours of the lenses, 
we expect the galaxy samples to be similar — but not identical 
— a nd the differences betw een the samples could have some ef- 
fect. iLeauthaud et al.l 1 120120 did not split their sample in red and 
blue, which is why we show the same constraints in both pan- 
els of Figures [12] and [T3] Further variations between the anal- 
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Figure 12. Comparison between four different datasets. The left (right) panels show the measured halo mass as a function of luminosity (stellar mass), and 
the top (bottom) panels show the results for red/early-type (blue/late-type) galaxies. The datasets used are all based on galaxy-galaxy lensing analyses with 
solid dots showing the CFHTLenS results from this paper. Also shown are halo masses measured using the RCS2 (open stars; VUl 1), the SDSS (open squares 
iMandelbau m et al. 2006) and COSMOS (solid band; Leauthaud et al. 2012). In the case of COSMOS we use the results from their lowest redshift bin. Also 
note that no distinction between red and blue lenses was made in the COSMOS analysis, so the same results are shown in both right panels. 



yses are discussed in more detail below. With these caveats in 
mind, we observe that all studies find similar general trends, with 
a halo mass that increases with increasing luminosity and/or stellar 
mass. It is also clear that blue/late-type galaxies tend to reside in 
haloes of lower mass than red/early-types do. The halo mass es- 
timates of blue galaxies presented in these studies are in excel- 
lent agreement. For the red galaxies, our mass estimates are con- 
sistent with those from VUllii and Mandelbaum et al.l ( 120061) ex- 
cept near Lr ~ lO^'^ h^o Lq, where they are 2-3cr lower. How- 
ever, as a function of stellar mass, our mass estimates of early-types 
broadly agree with th eirs. The halo masses of early-types also agree 
with the results from lLeauthaud et al.l (120121) at stellar masses be- 
low A/, ~ 10^^ hjQ Mq. At higher s tellar masses, the mass esti- 
mates are ~2ct lower than those from lLeauthaud et al.l ( |2012|). but 
we note that this is also the case for the halo masses from I VU 111 
and Mandelbaum et al. (2006). We will discuss this in more detail 
below. In general, a consistent picture of the relation between the 
baryonic properties of galaxies and their parent haloes is emerging 
from the four independent studies. 

Since our halo model is most closely related to that used by 

IVUllI (shown as open stars in Figures [T2l and [TJt, a detailed com- 
parison is more straight-forward compared to the other analyses. In 
IVUllL 1.7 X 10* lens galaxies were studied using the overlap be- 
tween the SDSS and the RCS2. The combination of the two surveys 
allowed for accurate baryonic property estimates using the spec- 
troscopic information from the SDSS, and a high source number 
density of 6.3 arcmin^^ owing to the greater depth and better ob- 
serving conditions of the RCS2 compared to the SDSS. Because 



we use photometric redshifts for our analysis our lens sample is 
more than sixty times that of VUll, reflecting the small fraction 
of galaxies that have spectroscopic redshifts determined by SDSS. 
The even greater depth of the CFHTLenS compared to the RCS2 
means that o ur sourc e density is a factor of 1.7 higher. Furthermore, 
in contrast to lVUlll we have individual redshift estimates available 
for all our sources. The increased number density and redshift res- 
olution in our analysis results in significantly tighter constraints on 
the relations between halo mass and luminosity, and between halo 
mass and stellar mass. 

As evide nced by Figure [T2| our halo masses agree well with 
those found bv lVUlll in general, though our halo mass relations are 
shallower; for red lenses we measure a power law slope for the rela- 
tion between halo mass and luminosity of LSGIoq^, ^^"1 between 
halo mass and stellar mass of 1.491qq4, while I VU 111 find slopes 
of 2.34lQig, and 1.9 ±0.1, respectively, using the same power 
law definitions. The general trend with stellar mass of a dec reasing 
baryon conversion efficiency for red lenses was observed bv lVUllI 
as well, but they were unable to discern a trend in their late-type 
sample. There are some differences between the analyses which 
should be noted, however. As mentioned above, we divide o ur lens 
sample in a red and blue one based on the SED type, while I VU 111 
use the brightness distribution profiles to separate their lenses in a 
bulge-dominated and a disk-dominated sample. Even though the 
resulting samples are expected to be fairly similar, they are not 
identical. As the mass-to-luminosity ratio of galaxies strongly de- 
pends on their colour, even small colour differences between the 
samples could result in different masses. This may explain why our 
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halo mass estimates o f the re d le nses at the high luminosit y end 
are lower than those of lVUllI and lMandelbaum et alj 1 120061) , who 
both use identical galaxy type separation criteria and whose masses 
agree in this regime. The difference is smaller for the stellar mass 
results, providing further support for this hypothesis. 

Furthermore, in our halo model we account for t he bary onic 
mass of each lens, something that was not done in IVUlll Re- 
moving the baryonic component from our model, we find that the 
masses for some bins increase by as much as 15%. It may appear 
counter-intuitive that including a baryonic component with a mass 
which is of order 5% of the total mass should result in such a sig- 
nificantly lowered halo mass estimate. The explanation lies in the 
halo model fitting, and specifically in the way the satellite fraction 
is allowed to vary. Adding a baryonic component on small scales 
will result in a lowered central halo mass. The central halo pro- 
file reaches further than the baryonic component however, and thus 
power on intermediate scales is also diminished. To compensate 
for this loss of power, the halo model will increase the satellite 1- 
halo term by increasing the satellite fraction, which also increases 
the stripped satellite halo term, lowering the central 1-halo term 
further until an equilibrium is reached. These mechanisms are il- 
lustrated for red galaxies in luminosity bin L4 in Figure [14] where 
we have allowed halo mass, satellite fraction and stellar mass frac- 
tion to vary simultaneously for both panels. This figure also makes 
clear the degeneracies introduced to the halo model if the stellar 
mass is left as a free parameter in addition to halo mass and satellite 
fraction. Higher luminosity or stellar mass bins are more severely 
affected by this effect than the lower end due to the lack of a promi- 
nent satellite 1-halo feature. To study the effect on the best-fit power 
law parameters, we re-fit the halo models excluding the baryonic 
component. The resulting slope and amplitude of the pow er law 
do not change significantly. Hence this does not explain whv lVUllI 
find a steeper slope than we do. We note, however, that our baryonic 
component only accounts for the stars in the lens and not for exam- 
ple the hot gas. The influence of feedback on the gas distribution in 
galaxies is a complicated issue which may also affect our re s ults, a s 
discussed in lvan Daalen et alj ( 120 Uh and lSemboloni et alj j2011r) , 
but it is an effect which we do not attempt to model here. However, 
as future lensing surveys grow more powerful and the data allows 
for greater accuracy this will be an important effect to study. 

Another factor to take into account is the fact that we limit 
our lens samples to redshifts of 0.2 ^ zicns ^ 0.4 keeping 
our me an lens redshift fairly stable at (zicns) ~ 0.3. This is not 
done in I VU 111 and as a result the median redshift of our lower 
lu minosit y or stellar mass bins is higher than for the same bins 
in IVUllL with the opposite being true for the higher bins. Re- 
cent numerical simulations indicate that the relation between stel- 
lar mass and halo mass wi l l evolve with redshift (see for example 
IConrov & Wechsled 120091: iMoster et all 12010). Lower-mass host 
galaxies (M* < 10^^ Mq) increase in stellar mass faster than their 
halo mass increases, i.e. for higher redshifts the halo mass is lower 
for the same stellar mass. The opposite trend holds for higher- 
mass host galaxies (Af, > 10^^ Mq). As a result, the relation 
between halo mass and stellar mass (or an indicator thereof, such 
as luminosity) steepens with increasing redshift. This means that 
for the lower-luminosity bins, w here ou r redshifts are higher, we 
may measure a steeper slope than lVUlll and vice- versa for higher- 
luminosity bins. The effect is likely small, however, because of the 
relatively small redshift ranges involved. 

Finally we note that the lenses in the sample studied bv lVUllI 
are rather massive and luminous as only galaxies with spectroscopy 
are used. Our lens sample includes many more low luminosity and 
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Figure 13. Comparison between four different datasets, showing the ra- 
tio of measured halo mass to stellar mass as a function of stellar mass. 
The top (bottom) panels show the results for red/early-type (blue/late-type) 
galaxies. The datasets used are all based on galaxy-galaxy lensing analy- 
ses with solid dots showing the CFHTLenS results from this paper Also 
shown are halo masses measured using the RCS2 (open stars; IVUllh . the 
SPSS (open squares M andelbaum et al. 2006) and COSMOS (solid band; 
iLeauthaud et al.ll2012l) . In the case of COSMOS we use the results from 
their lowest redshift bin. Also note that no distinction between red and blue 
lenses was made in the COSMOS analysis, so the same results are shown 
in both panels. 
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Figure 14. Dependence of halo model fitting parameters halo mass M200 
and satellite fraction a on stellar mass, with /sm the fraction of true mean 
stellar mass used in the halo model and contours showing the 67.8%, 95.4% 
and 99.7% confidence intervals. The left panel shows that including a bary- 
onic component in the model (i.e. setting /sm = 1) will result in a signifi- 
cantly lower best-fit halo mass than not doing so (/sm = 0), and the right 
panel shows that the reason for this is an increased satellite fraction. In our 
analysis we keep the stellar mass component fixed at /sm = 1- 
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low stellar mass objects, however. Hence the difference in slope 
may be partly due to the fact that we probe different regimes, and 
that the relation between baryonic observable and halo mass is not 
simply a power law but t urns upward at high lum inosities/stellar 
masses, as the results from lLeauthaud et al.l ( 120121) suggest. 

Having compared our analysis to that of VUl l ], we now turn 
our attention to the comparison with the Mandelbau m et al.l ( 120061) 
analysis of 3.5 x 10^ lenses in the SDSS DR7, shown as open 
squ ares in Figures [12] and [13] Their lens sample is, similarly to 
the IVUllI sample, also divided into early- and late-type galaxies 
based on their brightness profiles. To allow for a comparison be- 
tween our results and theirs we first have to consider the differ- 
ences in the luminosity definition. Man delbaum et alj ( l200q) use 
absolute magnitudes which are based on a K correction to a red- 
shift of z = 0.1 and a distance modulus calculated using h = 1.0. 
Furthermore, their luminosities are correcte d for p assive evolution 
by applying a factor 1.6(2: — 0.1). However. lVUl ll convert their lu- 
mino sities, which are similar to ours, using the JMandelbaum et al.l 
( I2OO6) definition and find that for low-luminosity low-redshift sam- 
ples the difference between the two definitions is negligible. The 
more luminous lenses reside at higher redshifts and for them the 
correction is found to be greater, most likely due to the differ- 
ence in the passive evolution corrections. Since our lenses are con- 
fined to relatively low redshifts, and since the main difference be- 
tween luminosity defin itions is the passive evolut ion factor, we can 
compare our results to JMandelbaum et al.l (I2OO6 ) without correct- 
ing the lumi nosities. Our halo mass d efinition is also different to 
that us ed bv [Mandelbaum et alj ( 120060 though. [Mandelbaum et al.l 
( [2OO6 1) define the mass within the radius where the density is 180 
times the mean background density while we set it to be 200 times 
the critical density. The correction factor stemming from the dif- 
ferent definitions amounts to ~ 30%. Having corrected for this, 
our re sults are then very similar to those from [Mandelbaum et al.[ 
( [20060 . but the same concerns of object selection and baryonic 
contribution discussed above apply here as well. The relation that 
[Mandelbaum et al.[ ((2006) find between halo mass and luminosity 
for red lenses is shallower than the one found by iVUll . as dis- 
cussed therein, and are therefore more in agreement with our re- 
sults. For the stellar mass relation, however, they find a steeper 
power law slope, though this resul t is mo stly driven by their highest 
stellar mass bin as pointed out bv [VUll[ 

Finally, [Leauthaud et al.[ (|2012|) perform a combined analy- 
sis of galaxy-galaxy lensing, galaxy clustering and galaxy number 
densities using data from the COSMOS survey, shown as a solid 
band in the right panels of Figure[T2land in Figure[T3] For our com- 
parison we select the results from their lowest redshift bin, since its 
redshift range of 0.22 < z < 0.48 is very similar to t he redshift 
range used here. Contrary to the other datasets, Leaut haud et al.[ 
ilOlT i did not separate their lens sample according to galaxy type. 
The results shown in the top panel of Figures [T2I and [T3l are there- 
fore identical to those shown in the bottom panel. Note that at 
high stellar masses, their sample is expected to be dominated by 
red galaxies, and at low stellar masses by blue galaxies, as these 
are generally more abundant in the respective regimes (see Table 
[3]l. Fo r stellar masses lower t han 10^^ h^g Mt, the agreement be- 
tween [Leauthaud_etalJ ( [2012h and the other galaxy-galaxy lensing 
results is excellent for both galaxy types. For higher stellar masses, 
however, Leauthaud et al. (2012) find higher halo masses than what 
has been observed in the lensing only analyses discussed above. 
Thi s may be explained if a larger fraction of the galaxies used in 
the [Leauthaud et al.[ ( [2012[) analysis reside in dense environments 
and can be associated with galaxy groups and clusters such that 



their halo masses correspond to the total mass of these structures . 
This theory is corroborated by Figure 10 of [Leauthaud et al.[([2012h 
which shows that for large stellar masses, the ratio of stellar mass 
to halo mass is ver y similar to that determined for a set of X-ray lu- 
minous clusters in [Hoekstral ( [20070 , indicating that we are entering 
the cluster regime. Furthermore, the sampling variance is not taken 
into account in the COSMOS error range. This is likely to affect 
the higher stellar mass bins more because the number of objects 
there is sparse. Additionally, the results fr om the COSMOS analy - 
sis of X-ray selected groups presented in [Leauthaud et al.[ ([20100 . 
which is centred on a redshift sim ilar to ours and also shown in Fig- 
ure 10 of [Leauthaud et al.l ( [2012[) as grey squares, agree better with 
our results for higher stellar masses. We note, however, that an- 
other possibility is that the high stellar mass end constraints from 
[Leauthaud et al.[ ( [20 120 may be driven mainly by the stellar mass 
function rather than by the lensing measurements. This, combined 
with the differences in the two halo model implementations, could 
also contribute to the observed discrepancy. 

A further subtlety discussed in Section |43| is that the satellite 
fraction of galaxies with high stellar masses is not well constrained 
by galaxy-galaxy lensing only. Since the satellite fraction and halo 
mass are weakly anti-correlated (see VUl 1), our halo masses may 
be slightly underestimated if the satellite fractions are too high. Fur- 
thermore, the modelling of the shear signal from satellites in this 
mass range is a bit uncertain as they may have been stripped of 
more than the 50% of their dark matter we have assumed so far, 
and this could also have some effect. However, we estimate that 
these modelling uncertainties only have a small effect on our best- 
fit halo masses, and that it is not sufficient to explain the differences 
between the results. 



7 CONCLUSION 

In this work we have used high-quality weak lensing data pro- 
duced by the CFHTLenS collaboration to place galaxy-galaxy lens- 
ing constraints on the relation between dark matter halo mass and 
the baryonic content of the lenses, quantified through luminosity 
and stellar mass estimates. The combination of large area and high 
source number density in this survey has made it possible to achieve 
tighter constraints compared to previous lensing surveys such as 
the SDSS, COSMOS or the RCS2. We also extended our study to 
lower stellar masses than have been studied before using a halo 
model such as the one described here. 

In this paper we have included a halo model constituent which 
was neglected by most earlier implementations: the baryonic com- 
ponent. Since the lensing signal is a response to the total mass of a 
system, it is essential to account for baryons in order to not overes- 
timate the mass contained in the dark matter halo. We have shown, 
however, that care has to be taken when including a baryonic com- 
ponent since doing so has a greater impact on the fitted halo mass 
than one might naively expect due to the complicated interplay be- 
tween stellar mass, satellite fraction and halo mass. 

As luminosity and stellar mass increases, the halo mass in- 
creases as well. For red lenses, the halo mass increases with 
greater baryonic content at a higher rate than for blue galaxies, 
independent of whether the measure of baryonic content is lu- 
minosity or stellar mass. The two measures thus produce com- 
parable results. For each we fit power law relations to quan- 
tify the rate of increase in halo mass. We find a best-fit slope 
of 1.56to°ot ™d a normalisation of 1.26lo;;;g x lO" /i^ M© 
for a fiducial luminosity of Lfld = 10^^ h^g Lq for red galax- 
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ies, while for blue galaxies we find a slope of 0.731q Qg 
normalisation 0.16 ± 0.03 x 10^^ h^Q Mq. The power law rela- 
tion between stellar mass and halo mass has a slope of IAQ'^^q^ 



and a noiTnalisation of l.SOlJJ'Qg x 10^^ h^^ Mq for a fiducial 



mass of Mfld = 2 x 10^^ h^g Mq for red galaxies, and for blue 
galaxies we find a slope of 0.831q'q4 and a normalisation of 
0.88 ± 0.13 X 10^^ /itoMq. 

For our blue galaxy selection, the satellite fraction is low 
across all luminosities and stellar masses considered here. The sig- 
nal at large scales for these samples is also generally low in the 
lowest luminosity and stellar mass bins, indicating that these galax- 
ies are relatively isolated and reside in less clustered environments 
than the red galaxies do and that we may be overestimating the 
galaxy bias for these samples. At low luminosity/stellar mass, a 
considerable fraction of red galaxies are satellites within a larger 
dark matter halo. This fraction decreases steadily with increasing 
luminosity or stellar mass. In general, the satellite fractions show 
that at these redshifts the galaxies in denser regions are mostly red 
while for the same luminosity or stellar mass isolated galaxies tend 
to be bluer and thus star forming. This indicates that the star for- 
mation history of galaxies differs depending on the density of the 
environment they are residing in. Furthermore, since blue galaxies 
have a fairly shallow relation between halo mass and baryonic mass 
there are a range of stellar masses possible for a given halo mass. 
This is not true for red galaxies, which therefore have a narrower 
range of baryon conversion efficiencies for a given halo mass than 
blue galaxies. 

Another finding of this work is that for faint and low stellar 
mass blue galaxies, the amplitude of the lensing signal at projected 
separations larger than ~2ftf(,^Mpc is lower than the corresponding 
best-fit halo model. For the red galaxies, the halo model fits the 
data well over all scales. This could indicate that while the bias 
description works well for red galaxies, it is not optimal for blue 
galaxies. If this is the case, then the environments in which the two 
samples reside are radically different and the difference will have to 
be taken into account in the future. Alternatively, the discrepancy 
could be caused by other choices that affect the 1-halo to 2-halo 
transition regime in our halo model implementation. Currently, we 
do not have enough data to favour or rule out either scenario, but 
we plan to explore this further in upcoming works. 

The relations between baryonic content indicators and dark 
matter halo mass presented in this work, as well as the dependence 
of the satellite fraction on luminosity and stellar mass, improve our 
understanding of the mechanisms behind galaxy formation since 
they provide constraints that can be directly compared to numerical 
simulations that model different galaxy format ion scenarios. With 
currently ongoing (for instance DESnor KiDSi lde Jong et alj2013r) 
and p lanned (such as LSSlQ HSQj or Eucliqj; iLaureiis et al.l 
l2011i) surveys, weak lensing analyses will become yet more power- 
ful than the one presented in this paper. In preparation for the future 
there are therefore several sources of uncertainty that should be in- 
vestigated. As mentioned above, the galaxy bias description may 
not be optimal for blue lenses and with future data this bias can 
likely be constrained directly using galaxy-galaxy lensing observa- 
tions. Recent simulations have also indicated that there is a redshift 
evolution of the halo mass relations, and this evolution can be stud- 
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led with weak lensing (see lChoi et"Zll2012l : [Hudson et alj|2013h . 
Other possible improvements to the halo model used here include 
studies of the distribution of satellites within a galaxy dark matter 
halo, a more accurate description of the regime where the 1-halo 
term and 2-halo term overlap (i.e. halo exclusion), and investiga- 
tions into the stripping of satellite haloes. The analysis presented 
in this paper is a significant advance on recent analyses, and with 
future surveys we will be able to use galaxy-galaxy lensing to study 
the connection between baryons and dark matter in exquisite detail. 
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APPENDIX A: CORRECTIONS FOR SIGNAL 
CONTAMINATION 

Al Photometric redshift bias correction 

Though the quality of the CFHTLenS photometric redshift esti- 
mates is high, there is still a small bias present due to the inherent 
limitations of template- based Bayesian methods, as discussed in 
iHildebrandt et al.l ( l2012h . This bias will affect not only the redshift 
itself, but also the derived quantities such as luminosity and stellar 
mass. Since our lenses reside at relatively low redshifts we there- 
fore have to coiTect our lens redshifts and derived quantities for 
this bias in order to achieve accurate object selection for our dark 
matter halo relations. Additionally, if this bias is not corrected for, 
the angular separations between lenses and sources will be altered, 
causing a coherent shift in the lensing signal radial binning. The re- 
sulting halo model fit will then also be affected, further illustrating 

the importance of this correction. 

Following IHildebrandt et alj ( |2012|) . we perform our correc- 
tion using spectroscopic redshifts in the overlap with the VI- 
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MPS VLT Deep Survey (VVDS; lLe Fevre et al.|2od5l:lGarilli et al 



200S ) , the DEEP2 galax y redshift surv ey dDavis et alj 



2007; Newman et al. 2012) and the SDSS ('Eisenstein et alj 



] 



2003 



2001 



Strauss et al. 2002). The completene ss of this spec troscopic sam- 
ple is s hown in lLe Fevre et al.l l l2005l . Figure 16) and lNewman et al.l 
1 120 I2I Figure 31). To ensure a completeness of at least 80%, we 
select only lenses with magnitude i^g < 23, as mentioned in Sec- 
tion l2.2l Since the bias is a function of magnitude and galaxy type, 
we start by splitting our sample in red and blue subsamples via 
their photometric type (as described in Section l2Ti and use several 
magnitude bins. We then quantify the bias in each bin by fitting a 
straight line of the form 



(^spoc - 0.3) = a(zphot — 0.3) + b 



(Al) 



where Zspcc is the spectroscopic redshift from 
VVDS/DEEP2/SDSS, Zphot is the CFHTLenS photometric 
redshift estimate, a is the slope and b is the offset. The pivot point 
of 0.3 roughly corresponds to the mean redshift of our lens sample, 
though the correction is insensitive to this number. The slope a 
is fit simultaneously in all magnitude bins but allowed different 
values for red and blue samples, while the offset b is allowed to 
vary between both type and magnitude bins. Keeping the slope 
fixed allows for a more robust estimate for the bias, though we have 
verified that allowing it to vary has negligible impact on the results 
in practice. The resulting fit parameters are shown in Table lAll 
Note that there is no correction performed for red galaxies beyond 
a magnitude of Jab = 22 since we do not use fainter red lenses 
(see Section [2!2l l. 

We then use these fit parameters to correct our lens photomet- 
ric redshift estimates in the range 0.2 ^ 2; ^ 0.4. Calculating the 
luminosity distances and estimating the k-corrections correspond- 
ing to the original and corrected redshifts using the g' — r' colours 
of the galaxies, we adjust the absolute magnitudes accordingly. We 
further derive new stellar mass estimates by scaling them to their 
new luminosities assuming a constant (pre-correction) stellar mass- 
to-luminosity ratio. The impact on the red galaxy properties is neg- 
ligible, but for blue galaxies the correction is larger with the aver- 
age luminosity and stellar mass increasing by ~ 12%. We therefore 
proceed to use the corrected quantities in our luminosity and stellar 
mass analyses (see Sections|4]and|5]l- 

The sources will also be affected by photometric redshift bias, 
but its impact on the measured halo masses is expected to be much 
smaller than the effect of the lens redshift bias. To confirm this hy- 
pothesis we shift all sources by a constant bias of 2% and redo the 
analysis of Sections[4]and [5] This bias value is mos t likely slightly 
larger than necessary (see iHildebrandt et al.ll2012l Figure 4), but 
the resulting halo masses agree with the original halo masses within 
1(7. We therefore do not need to correct our sources for photometric 
redshift bias. 



A2 Photometric redshift scatter correction 

Before interpreting the luminosity results we have to take into ac- 
count the effect of Eddington bias JEddingtonll 1 9 1 31) . The precision 
of our photometric reds hifts is high with a scatte r of az ~ 0.04 for 
both lenses and sources JHildebrandt et al.l2012h . but never the less 
the errors on the redshift estimates have to be taken into account. If 
the true redshift differs from the estimated one, this will affect all 
derived quantities. An underestimated redshift, for example, would 
cause the estimated absolute magnitude to be fainter than the true 
absolute magnitude and the lens would be placed in the wrong lu- 
minosity bin. As can be seen in Figure|4]there are more faint objects 




Figure Al. Con'ection factor as a function of luminosity induced through 
inaccuracies in the photometric redshift estimates. The dark purple solid 
(light green dotted) line with dots (triangles) shows the scatter correction 
factor for the red (blue) lens sample. The eiTor bars show the scatter between 
10 lens catalogue realisations. 



Table A2. Photometric redshift scatter correction factors applied to ob- 
served halo masses in each luminosity bin (see Section O for red and blue 
lenses. These factors connect both for scatter due to redshift en'ors, and for 
the fact that the observed halo mass does not necessarily correspond to the 
mean halo mass. 
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than bright, which means that more objects will scatter from fainter 
bins into brighter bins than the other way around. This will lower 
the lensing signal in each bin and bias the observed halo mass low, 
and the amount of bias will be luminosity dependent. To estimate 
the impact of redshift scatter we create a simulated version of the 
CFHTLenS as follows. We fit an initial power law mass-luminosity 
relation of the form (see Equation[TO] Section HTTt 



M200 = Mo,. 



L 

ifld 



(A2) 



10 /I7Q ir',0- 



to the raw estimated halo masses, with Lad 
We then use this relation to assign halo masses to our lenses. Split- 
ting the resulting lens catalogue in the usual magnitude bins for 
the red and blue sample separately we obtain our 'true' halo mass 
for each bin. Constructing NEW haloes from these halo masses at 
the photometric redshift of the lenses, we create mock source cata- 
logues with the observed source redshift distribution but with sim- 
ulated shear estimates with strengths corresponding to those which 
would be induced by our lens haloes. We then scatter the lenses 
and sources using the full redshift probability density function, split 
the lens catalogue according to the scattered magnitudes and mea- 
sure the resulting signal within 200 h^^ kpc of the lenses using 
our scattered shear catalogue. We only use the small scales for our 
mass estimate to avoid complications due to insufficient treatment 
of clustering since on these scales only the central 1-halo signal 
is relevant, and we force our satellite fraction to zero to obtain a 
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pure NFW fit. This way we obtain the 'observed' halo mass for 
each magnitude bin. The 'observed' halo mass is then compared to 
the 'true' value for each bin. To increase the statistical precision of 
the correction we determine the average of 10 lens catalogue real- 
isations. Since the starting point is a perfect signal, the number of 
realisations given the area is adequate to retrieve the correction fac- 
tor. This correction simultaneously accounts for all the effects re- 
sulting from any photometric redshift scatter in our analysis, such 
as the scattering of lenses between luminosity bins, the effect on 
the lens and source redshift distributions, the smoothing of the sig- 
nal due to mixing of the projected lens-source separations, and the 
non-linear dependence of the critical surface density Ecrit on the 
lens and source redshifts. Note that the errors on the correction fac- 
tors indicate only the propagated photometric redshift uncertainty, 
and even though they are small compared to the errors on the shear 
measurements, we have included them in our final error budget. 
The error on the correction factor does not include the uncertain- 
ties of the input parameters. However, we expect these additional 
uncertainties to be negligible compared to the errors on the halo 
masses (see the discussion in Appendix lA3t . 

The results from this test are shown in Figure IaTI The quality 
of our photometric redshifts is high which means that the correc- 
tion factor is small overall, reaching only ~ 30% for a luminos- 
ity of Lj./ r^ 2.5 X 10^^ h^Q Lq. Here the contamination is largest 
due to the shape of the luminosity function causing a larger frac- 
tion of low luminosity objects to scatter into the higher-luminosity 
bin. For our faintest red luminosity bin the correction is also 
~ 30%, in this case caused by larger errors in the photometric red- 
shift estimates. The correction factor is less than unity for lower- 
luminosity bins due to the turn-over of the distribution of red lenses 
at M^i ~ —21.2 (see Figure^. The small correction factor for blue 
lenses is due to their flatter mass-luminosity relation (see Figure|6f. 
Because of the relative insensitivity of halo mass to changes in lu- 
minosity, minor errors in luminosity measurements due to photo- 
metric redshift inaccuracies will not strongly affect the halo mass 
estimate. The process described in this appendix could in princi- 
ple be iterated over, starting from the fittin g of a mass-lum i nosity 
relation, until convergence is reached. Since iHoekstra et alj j2005r) 
find that different choices for that relation yield similar curves, we 
choose not to iterate further. 

We also have to corre ct our l uminosity bins for a second scat- 
ter effect. As discussed in I VU 111 (Appendix B), the observed halo 
mass does not necessarily correspond to the mean halo mass in a 
given bin since the halo masses in that bin are not evenly distributed 
and the NFW profiles do not depend linearly on halo mass. The 
distribution within each bin generally follows a log-normal distri- 
bution, and to correct for this effect we follow a similar procedure 
as the one outlined in Appendix I A3 1 with the difference that we do 
not scatter the luminosities as we have already corrected for that by 
accounting for the error in photometric redshift. We stress that this 
is an intrinsic effect unrelated to any measurement errors. The full 
correction factor, taking into account both scatter eifects discussed 
here, is shown in Table lA2l 

The general procedure outlined in this appendix is repeated 
for the stellar mass bins, though we use the stellar mass-halo mass 
relation to assign halo masses to the mock lens catalogue, and then 
bin the lenses according to stellar mass rather than luminosity. In 
this case we do not use the resulting correction factor, but we do 
include the errors on said correction factor in our error budget in 
order to account for the above-mentioned effects in our stellar mass 
results. The correction factor itself, however, only conveys the im- 
pact of photometric redshift uncertainties, and not the additional 
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Figure A2. Mass correction factor as a function of stellai' mass induced 
through inaccuracies in the stellar mass estimates. The dark purple solid 
dots (light green open triangles) show the correction factor for the red (blue) 
lens sample. As discussed in the text, the dot-dashed lines show the correc- 
tion factors if stellai' mass eiTors of 0.2 dex are assumed, rather than the 
default 0.3 dex. and the dotted lines show the correction factors derived 
using stellar mass errors of 0.4 dex. 



Table A3. Bin scatter correction factors applied to observed halo masses in 
each stellar mass bin (see Section|5j for red and blue lenses. These factors 
correct both for scatter due to stellar mass en'ors, and for the fact that the 
observed halo mass does not necessarily correspond to the mean halo mass. 
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effects influencing the stellar mass errors. The scatter due to stellar 
mass errors is accounted for following the method described in the 
next appendix, and applying this correction factor as well would 
therefore amount to double-counting. 

A3 Stellar mass bin scatter correction 

In a process similar to the scatter in luminosity, objects will scatter 
between stellar mass bins due to errors on the stellar mass esti- 
mates. Though objects scatter randomly according to their individ- 
ual stellar mass eiTors, the net effect will be to scatter lenses from 
greater abundance to lower according to the stellar mass function 
(SMF). Because the stellar mass function declines steeply at higher 
stellar mass bins, these will be more severely affected by low-mass 
object contamination. As a result the observed lensing mass in the 
highest stellar mass bins will be biased low (see Appendix A in 
IVUID . Additionally, the lensing halo mass estimates will be af- 
fected by the fact that the observed halo mass does not necessarily 
correspond to the mean halo mass in a given bin, as discussed in 
Appendix I A2 1 

To assess the impact of both these effects simultaneously we 
follow a procedure similar to the one used to correct for redshift 
scatter, as described in the previous appendix. We start by fitting 
an initial power law halo mass-stellar mass relation using the raw 
observed lensing halo mass. Drawing a large number of simulated 
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lens galaxies from the stellar mass function, we take these stellar 
masses to be the true unscattered values and assign a halo mass ac- 
cording to the fitted halo mass-stellar mass relation. As described 
above, this halo mass will be distributed w ithin th e stellar mass bin 
according to some distribution. Following I VU 1 ll we therefore cor- 
rect the halo mass for this effect by drawing from a log-normal dis- 
tribution with a mean given by th e original halo mass and a width 
determined by Mor e et al.l 1 1201 11) . We now know the 'true' mean 
halo mass for each bin. Using the resulting simulated lens catalogue 
we create a source catalogue with shears determined analytically. 
We then scatter the lenses assuming a Gaussian error distribution 
with a width of 0.3 dex as appropriate for our stellar mass errors 
(see Section im to create a new simulated lens catalogue, this time 
containing 'observed' stellar masses. Dividing this 'observed' lens 
catalogue according to the usual stellar mass bins for red and blue 
samples separately, we measure the signal in the simulated shear 
catalogues and again fit an NFW profile. This way we obtain the 
'observed' halo mass for each stellar mass bin. By taking the ratio 
of simulated 'observed' to 'true' halo mass we aiTive at the cor- 
rection factor for stellar mass scatter as shown in Figure IA2I We 
can now apply this factor, as quoted in Table I A3] to our halo mass 
estimates to correct for the scatter between stellar mass bins, and 
for the fact that the observed halo mass does not correspond to the 
mean halo mass, simultaneously. 

The correction factor is relatively sensitive to the adopted 
value of the stellar mass error, particularly in the regime where 
the stellar mass function is steep. Therefore, in addition to the 
correction factor used, we also show in Figure lA2l the correction 
factors obtained if we adopt a stellar mass error of 0.2 dex or 
0.4 dex instead, covering the plausible range of values that the 
stellar mass error could take. This illustrates how the correction 
factor coherently shifts if the stellar mass error is different from 
what we assume. For S8 of the red lenses, the change is largest, 
with an increase (decrease) of the coiTection factor by ~ 50% 
for 0.4 dex (0.2 dex), respectively. We do not use the plausible 
range of correction factors as the error on the correction, since a 
different stellar mass error would only lead to a coherent shift of 
all the correction factors and hence of the corrected halo masses. 
This property of the correction factors would be lost, and the error 
bars on the halo masses would be severely overestimated, caus- 
ing an unjustified loss of information. However, for complete- 
ness, we note that the best fit power law normalisation and slope 
are 1.00 ± 0.05 x lO" ft^ A/© (Omto\'l x lO" /i^ ^^0) and 
1.35 ± 0.05 (0.78t'oM) for red (blue) lenses when we adopt 
a stellar mass error of 0.2 dex, and l.eslg n x lO" ^^m -^0 
(1.17 ±0.18 X 10"/ito M0) and IMto.oilO-STllo^ol) for red 
(blue) lenses for a stellar mass error of 0.4 dex. 

Additionally, the correction factor has some error due to the 
uncertainties of the other input parameters, such as in the adopted 
power law relations, the stellar mass functions, and the scatter in 
halo mass. IVUllI found that the correction is fairly insensitive to 
changes in the power law relation; using the power law obtained 
after the stellar mass scatter correction only changed the correc- 
tion factor by at most 4%. The impact here will be even smaller 
as the power laws are less steep, and we therefore ignore their ef- 
fect. Next, the stellar mass function is not the intrinsic stellar mass 
function as objects have already scattered. However, we cannot reli- 
ably obtain the intrinsic stellar mass function where it matters most, 
i.e. at the high stellar mass range, as the number of galaxies is too 
low. We therefore do not attempt to obtain the intrinsic stellar mass 
function, but rather note this as a caveat. Finally, we note that the 
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Figure Bl. Comparison of three data sets: the shear catalogues from 
~ 22 deg2 CFHTLS (pink open squares), the results from RCS2 (light 
green open stars) and our results (dark purple solid dots). The curves show 
the best-fit singular isothermal sphere for each dataset (with light green and 
pink nearly identical), and the grey triangles show the cross-shear from our 
results which should be zero in the absence of systematic errors. 



correction factor is insensitive to the adopted width of the halo mass 
distribution. 



APPENDIX B: TESTS FOR GALAXY-GALAXY LENSING 
SYSTEMATICS 

Bl Initial consistency analysis of the CFHTLenS catalogue 

In this study we use lenses and sources from the full 154 deg^ 
CFHTLenS catalogue. The accuracy of the CFHTLenS shears has 
been verified through several rigorous tests aimed at the study of 
cosmic shear (Heymans et al. 2012; Miller et al. 2013), but it is in- 
teresting to compare the galaxy-galaxy lensing signal with the re- 
sults from two previous analyses of a similar nature. The first is 
th e galaxy-galaxy lens ing analysis in the CFHTLS-Wide conducted 
bv Parker e t al. (2007), and the seco nd is based on the shear cata- 
logue from iVUll (see Section|6}. In Parker et al.l ( l2007h an area of 
~ 22 deg^ in i' was analysed, corresponding to about 14% of our 
area. Since they only had data from one band their analysis also 
lacked redshift estimates for lenses and sources, but they separated 
lenses from sources using magnitude cuts. The shear estimates for 
their sources were obtained using a versio n of th e technique in- 
troduced bv Kaiser. Squires. & Broadhurstl ( Il995n as outlined in 
iHoekstra et alj ( Il998h . These shear estimates were measured on a 
stacked ima ge rather than obt ained by fitting all exposures simulta- 
neously (see lMiller et alj|2013l for a discussio n on this). To avoid 
the strong PSF effects at the chip boundaries, IParker et al.l (120070 
limited their analysis to the unique chip overlaps. In contrast we 
are able to use all the data we have at our disposal. The data from 
lyUll is the subset of ~400 square degrees of the RCS2 with i'- 
band coverage, which is shallower than the CFHTLS and for which 
also no redshifts were available for the sources at the time of this 
analysis. 

To compare and contrast our lensing signal with these previ- 
ous works we mimic the analysis presented in Parker et al. (2007|) 



as closely as p ossible and apply th e same i -band magnitude cuts 
as employed in lParker et al.l ( l2007h . with 19.0 < Jar < 22.0 for 
lenses and 22.5 < i'ab < 24.5 for sources. IParkeretaJl (|2007|) 
boost their signal to correct for contamination by sources that are 
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physically associated with the lens, and we apply the same coiTec- 
tion factor to our values. The resulting galaxy-galaxy signal, scaled 
with the angular diameter distance ratio (/3) — (Dis/Da) = 0.49 
from Parker et al. ( 200 7), is shown as dark purple solid dots in Fig- 
ure lBll We also re-ana lysed the original shear catalogues used for 
the lParker et al.l ( l2007h analysis with the res ults sho wn as pink open 
squares in Figure IB II The signal from the IVUllI shape measure- 
ment catalogues of the RCS2 is obtained using a source selection 
of 22 < r' < 24 instead because the limiting depth in i' is 23.8 
for the RCS2. The measurements are also corrected for contamina- 
tion by physically associated sources, as described in V UllL and 
scaled with (/?) — 0.30 which is determined by integrating over 
the lens and source redshift distribut ions that were o btained from 
the CFHTLS "Deep Survey" fields l lllbert et al.ll200 6). The mea- 
surements are shown as light green open stars. Figure IbTI shows 
that the lensing signals generally agree well. We fit an SIS profile 
to the shear measurements that have been scaled by (/3) on scales 
between 7 and 120 arcsec, and find a scaled Einstein radius of fi; = 
0.277 ± 0.006 arcsec for our results, fS = 0.267 ± 0.011 arcsec 
for the lParker 61^11 ilOOlh measurements and rg = 0.262 ± 0.007 
arcsec for lVUl iL which are broadly consistent. 

The best-fit SIS profile corresponds to a velocity dispersion 
of ay = 97.9 ± 1.0k ms~^, which is low er than the ay = 
132 ±10 k m s~^ quoted in Parker et al.l ( l2007h . However, using the 
re-analvsed lParker et al.l (l2007b shear catalogue we find a velocity 
dispersion of cr„ — 96.6 ± 2.0kms~^. For the IVU 111 results we 
find a velocity dispersion of Ui, = 95. 4±1.3kms~^, slightly lower 
but in reasonable agreement with our results. Note that there are 
various small differences between the analyses, such as different 
effective source redshift distributions and different weights applied 
to the shears. Additionally we use the multiplicative bias correction 
factor for our measurements, while the other works did not have 
such a correction. All these differences could have small but non- 
negligible effects on the results. The dis crepancy with the velocity 
dispersion quoted in jParker et al.l ( |2007|) remains unexplained, but 
we conclude that the shear estimates are in fact fully consistent. 



B2 Seeing test 

iMiller et al.l 00131) isolated a general multiplicative calibration fac- 
tor as a function of the signal-to-noise ratio and size of the source 
galaxy, m{psN,r), using simulations. To confirm the successful 
calibration of the CFHTLenS shears in the context of galaxy-galaxy 
lensing, we study how a shear bias relates to image quality. In gen- 
eral a round PSF causes circularisation of source images which in 
turn can cause a multiplicative bias of the measured shapes if it is 
not properly corrected for Such a systematic would depend on the 
size of the PSF. Assuming that the systematic offset due to PSF 
anisotropy is negligible (a fair assumption given our correction for 
spurious signal around random lenses; see Sect ion |3JJand the de- 
tailed analysis of PSF residual errors in lHevman s et al. (2012)), and 
assuming that the shapes of very well resolved galaxies can be ac- 
curately recovered, the observed average shear in a galaxy- galaxy 
lensing azimuthal distance bin is related to the true average shear 
via 
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where 7°'^'^ is the observed shear, 7*'"° is the true shear, r, is the 
PSF size, ro is the intrinsic (Gaussian) size of the galaxy and A4 
is a value close to zero representing the multiplicative bias. The 
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Figure B2. Galaxy-galaxy lensing signal quantified thi'ough the best-fit Ein- 
stein radius (see Equation lB2) as measured in each of 6 seeing bins, accord- 
ing to Table lBl] The solid line shows the best-fit model using Eg nation IB 3 1 
while the dashed line shows the average Einstein radius assuming no bias. 



Table Bl. Details of the seeing bins. 



Sample 


-'Vflclds 


(r*) [arcsec] 


6^ [arcsec] 


<T6)e 


PI 


27 


0.50 


0.053 


0.005 


P2 


23 


0.57 


0.044 


0.006 


P3 


33 


0.62 


0.050 


0.005 


P4 


38 


0.67 


0.047 


0.005 


P5 


28 


0.72 


0.040 


0.006 


P6 


36 


0.80 


0.049 


0.005 



particular dependence on PSF size is the result of a ful l moments 
analysis (see for example lPaulin-Henriksson et al.l2008l) . 

Since the bias depends on the size of the PSF relative to the 
size of galaxies, data with a spread in seeing should enable us to 
determine the bias M directly from the data, thus allowing us to 
deduce the true performance of the shape measurement pipeline. 
The CFHTLS images have such a spread, with the best seeing being 
0.44 arcsec and the worst being 0.94 arcsec, and therefore provides 
us with a neat way of determining this bias. Since at small projected 
separations from the lens, the tangential shear signal is generally 
well described by an SIS profile: 
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where 6 is the distance to the lens and 6e is the Einstein radius, we 
therefore have a simple relationship between the observed Einstein 
radius and the true one: 
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By measuring the Einstein radius of the average lens as a function 
of seeing we can therefore determine both the true Einstein radius 
and the performance of the shape measurement pipeline. 

We select our lenses in magnitude and redshift as described 
in the main paper (Section [2. 2t . though we do not distinguish be- 
tween red and blue galaxies, and we also split our data according 
to Table IbTI Dividing the data according to image quality in this 
way may imply some minor selection effects, such as redshift and 
magnitude estimates being less accurate for worse seeing and thus 
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PSF. Since great care ha s been taken to correct for such effects (see 
iHildebrandt et al.ll2012l) we will assume here that the lens samples 
are comparable between seeing bins. Having selected our lenses, 
we measure the galaxy-galaxy lensing signal in each seeing bin 
and fit an SIS to the innermost 200 h^^ kpc. By fitting only small 
scales we avoid the influence of neighbouring haloes. The results 
are shown in Figure lB2l and quoted in Table IB ll We then fit the re- 
lation described by Eauation lB3l to the resulting Einstein radii and 
find a value of A^ = —0.071 ± 0.075. This is consistent with no 
bias, a fact which is further illustrated in Figure lB2] the data points 
agree with an average Einstein radius of 0.058 ± 0.003, shown as 
a dashed line. 



APPENDIX C: DETAILED LUMINOSITY BINS 

In this Appendix we show the decomposition of the best-fit halo 
model for red (Figure ICTt and blue (Figure IC2t lenses, split in lu- 
minosity according to Table [T] Showing the full decomposition is 
highly informative because it highlights some of the major trends 
and clarifies which effects dominate in each case. 

The baryonic component based on the mean stellar mass in 
each bin (dark purple dot-dashed line) becomes more dominant for 
higher luminosities, but the luminous size of the lenses also in- 
creases, making measurement of background source shapes in the 
innermost distance bins difficult. Thus it is not possible to reli- 
ably constrain the baryonic component with our data. Yet the ef- 
fect of including the baryons in our model is an overall lowering of 
the dark matter halo profile (dark purple dashed) compared to the 
model without baryons. For the red lenses we see that a consider- 
able fraction of the sample at lower luminosities necessarily con- 
sists of satellite galaxies, since there is a clear bump in the signal 
at inteiTnediate scales which has to be accounted for. This satellite 
fraction continuously drops as luminosity increases, and simultane- 
ously becomes more difficult to constrain since the combination of 
the stripped satellite profile (light green dash-dotted) and satellite 
I -halo terms (light green dashed) becomes almost indistinguish- 
able from a single NEW profil e for hi gh halo masses. This effect 
was discussed in more detail in lVUlll Appendix C. 

For the blue lenses, the signal becomes very noisy for the two 
highest-luminosity bins due to a lack of lenses. These two bins are 
therefore discarded from the full analysis in Section |4] In general, 
blue galaxies produce a noisier signal than red galaxies for the same 
luminosity cuts. This could be because blue lenses are in general 
less massive, and there are fewer of them which results in a weaker 
signal and a lower signal-to-noise for most bins. It could also be an 
indicator that the physical correlation between stellar mass and halo 
mass is noisier for these lenses. We also notice that nearly all blue 
lenses are galaxies located at the centre of their halo, rather than 
being satellites. This is consistent with previous findings. It is pos- 
sible that satellite galaxies in general are redder because they have 
been stripped of their gas and thus have had their star formation 
quenched. It could also mean that most blue galaxies in our analy- 
sis are isolated; we have made no distinction between field galax- 
ies and galaxies in a more clustered environment. If blue galaxies 
are more isolated than red ones then the contribution from nearby 
haloes (dotted lines) would also be less. It is clear from Figure IC2l 
that the large scales are not optimally fit by our model, and isola- 
tion may be one of the reasons since we assume the same mass-bias 
relation for blue galaxies as for red. With current data it is not pos- 
sible to constrain the bias as a free parameter, but with future wider 
surveys this could be done. 
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Figure CI. Galaxy-galaxy lensing signal around red lenses which have been split into luminosity bins according to Table [T] and modelled using the halo 
model described in Section lX2l The black dots are the measured diiferential surface density, AE, and the black line is the best-fit halo model with the separate 
components displayed using the same convention as in Figure [3] Grey triangles represent negative points that are included unaltered in the model fitting 
procedure, but that have here been moved up to positive values as a reference. The dotted eiTor bars are the unaltered error bars belonging to the negative 
points. The grey squares represent distance bins containing no objects. 
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Figure C2. Galaxy-galaxy lensing signal around blue lenses which have been split into luminosity bins according to Table [T] and modelled using the halo 
model described in Section lX2l The black dots are the measured diiferential surface density, AE, and the black line is the best-fit halo model with the separate 
components displayed using the same convention as in Figure [3] Grey triangles represent negative points that are included unaltered in the model fitting 
procedure, but that have here been moved up to positive values as a reference. The dotted eiTor bars are the unaltered error bars belonging to the negative 
points. The grey squares represent distance bins containing no objects. 
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APPENDIX D: DETAILED STELLAR MASS BINS 

The decomposition of the best-fit halo model for red and blue 
lenses, divided using stellar mass as detailed in Table [5] is shown 
in Figures lDn and lD2l respectively. 

By construction the baryonic component amplitude (dark pur- 
ple dash-dotted line) increases with increasing bin number, and so 
does the dark matter halo mass (dashed lines). Note that with our 
stellar mass selections we push to smaller and fainter objects, so 
the objects in the three lowest-mass bins are on average less mas- 
sive and less luminous than the galaxies in the faintest luminos- 
ity bin. In these bins, nearly all red galaxies are satellites, while 
for higher stellar mass bins the satellite fraction diminishes, a be- 
haviour which is consistent with the trends we saw for luminosity 
(Appendix [Ct. For the higher stellar mass bins, as for the higher 
luminosity bins, the sum of the satellite stripped and 1-halo terms 
result in a profile which resembles a single NEW profile, making 
the satellite fraction more difficult to determine. For the blue lenses 
we run into the same issues for the highest mass bin as for the high- 
est luminosity bins; the number of lenses is too small to constrain 
the halo model and so the bin has to be discarded. Furthermore, 
the satellite fraction is low across all blue lens bins indicating that 
these lenses are most likely isolated, which is consistent with the 
low large-scale signal and with our findings for luminosity. 
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Figure Dl. Galaxy-galaxy lensing signal around red lenses which have been split into stellar mass bins according to Table [5] and modelled using the halo 
model described in Section lX2l The black dots are the measured differential surface density, AE, and the black line is the best-fit halo model with the separate 
components displayed using the same convention as in Figure [3] Grey triangles represent negative points that are included unaltered in the model fitting 
procedure, but that have here been moved up to positive values as a reference. The dotted eiTor bars are the unaltered error bars belonging to the negative 
points. The grey squares represent distance bins containing no objects. 
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Figure D2. Galaxy-galaxy lensing signal around blue lenses which have been split into stellar mass bins according to Table [3] and modelled using the halo 
model described in Section [3^ The black dots are the measured differential surface density, AE, and the black line is the best-fit halo model with the separate 
components displayed using the same convention as in Figure |3] Grey triangles represent negative points that are included unaltered in the model fitting 
procedure, but that have here been moved up to positive values as a reference. The dotted error bars are the unaltered error bars belonging to the negative 
points. The grey squares represent distance bins containing no objects. 



